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1 . 0  INTRODUCTION  AND  LITERATURE  SURVEY 


1.1  INTRODUCTORY  COMMENTS 

The  simulation  of  boundary  layer  transition  in  test  facilities 
would  be  a  useful  capability  to  have  in  the  development  of  higher 
performance  and  more  efficient  vehicles.  An  obstacle  to  our  developing 
the  capability  to  simulate  transition  has  been  the  freestream 
disturbances  in  the  test  facilities  which  are  not  completely  similar  to 
the  disturbances  in  the  atmosphere.  '  Another  obstacle  has  been  that  the 
surface  roughness  and  other  fine  geometric  details  of  the  subscaie 
models  are  not  completely  similar  to  those  details  on  the  full-scale 
flight  vehicle.  A  related  problem  is  our  inability  to  systematically 
develop  boundary  layer  trips  to  fix  transition  in  a  specified  region. 
However,  (dimensional)  duplication  or  (nondimensional)  complete 
similarity  of  the  disturbances  and  roughness  are  excessively  severe  and 
impractical  requirements.  The  practical  need  exists  to  identify  the 
significant  properties  of  the  freestream  disturbances  and  roughness 
which  influence  stability  and  transition. 

An  obstacle  to  identifying  those  properties  has  been  the  inability 
of  theory  to  link  freestream  disturbances  to  tne  instability  waves  and 
other  oscillations  in  an  initial-value  problem.  Another  obstacle  has 
been  the  inability  to  show  how  roughness,  either  isolated  elements  or 
distributed  roughness,  affects  the  initial  amplitudes  of  the  various 
oscillations.  When  analyzing  such  a  problem,  the  theoretician  must 
know  all  of  the  possible  solutions  so  that  the  energy  from  the 
freestream  disturbances  is  channeled  to  the  proper  waves. 

When  the  author  and  Dr.  Shunichi  Tsuge  studied  an  initial-value 
problem  in  space,  their  objective  was  to  describe  the  evolution  o± 
disturbances  in  the  downstream  direction  in  terms  of  a  superposition  ot 
solutions  of  the  Orr-Sommerfeld  equation.  The  set  of  stability  waves 
and  two  solutions  representing  vortical  freestream  disturbances  appear 
in  the  analysis.  However,  two  additional  mathematical  poles  appear  in 
the  analysis.  These  poles  indicate  that  two  additional  solutions  exist 
of  the  Orr-Sommerfeld  equation.  These  standing  wave  solutions  are 
generalizations  to  oscillations  previously  analyzed  by  Rogler  and 
Reshotko  (Refs.  1,2)'.  This  report  documents  tne  structure  of  tnese 
standing  waves.  The  initial-value  problem,  which  provides  tne 
framework  to  link  the  initial  conditions  with  tne  amplitudes  and  phases 
of  the  various  waves,  will  be  summarized  in  later  writings. 

Waves  which  do  not  travel  or  propagate  are  called  standing  waves . 
They  may  be  steady  or  unsteady.  Standing  waves  of  a  simplier  form  than 
considered  herein  have  been  studied  using  the  unsteady  boundary  layer 
equations  with  the  oscillating  freestream 


This  freestream  oscillation  does  not  depend  on  position.  Analytical 
solutions,  and  related  numerical  calculations  and  experiments  have  been 
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extensively  reported  in  the  literature,  as  surveyed  by  Rogler  and 
Reshotko  (Ref.  3)  and  Telionis  (Ref.  4).  Loehrke,  Morkovin,  and  Fejer 
(Ref.  5)  reviewed  the  stability  and  transition  literature  in 
oscillating  boundary  layers.  The  author  believes  that  this  simple 
standing  wave  is  a  limiting  case  of  a  more  general  standing  wave  with 
far  broader  applicability. 

Before  we  examine  the  more  general  case,  however,  we  wish  to 
illustrate  how  an  oscillating  freestream  (1.1)  interacts  with  a 
parallel-flow  boundary  layer.  Since  the  x-derivatives  vanisn,  tne 
momentum  equations  for  small-amplitude  fluctuations  reduce  to 

-  o  >  =  -G*  "Aese -%/f>  =■  -  ra>fae‘a,t  ^  ^ 

The  reduced  form  is  merely  the  one-dimensional,  unsteady  diffusion 
equation  with  a  known  forcing  term,  px/p.  The  longitudinal  velocity 
fluctuation,  u,  is  uncoupled  from  both  the  mean  flow  and  tne 
fluctuation  velocity  normal  to  the  plate.  It  is  merely  the  flow  driven 
by  the  exterior  pressure  gradient  (which  oscillates  in  time)  and 
altered  by  the  diffusion  of  vorticity  to  and  from  the  wall  so  that  tne 
no-slip  condition  is  satisfied. 

However,  another  form  of  standing  waves  is  possible  in  the  form  of 
waves  which  grow  or  decay  exponentially  in  space,  as  well  as 
oscillating  sinusoidally  in  time.  Generally  these  waves  have  both  u 
and  v  velocity  fluctuations  and  are  coupled  with  the  mean  boundary 
layer  profile.  To  our  knowledge,  these  waves  have  not  be  studied  in 
any  shearing  layers  before,  although  they  have  appeared  in  several 
studies  with  uniform  mean  flows  which  are  summarized  in  the  following 
sections.  These  oscillations  are  believed  to  arise  in  many  practical 
situations,  and  have  a  role  in  the  initial-value  problem  in  space  for 
disturbances  described  by  the  Orr-Sommerfeld  equation.  These  waves 
also  have  a  role  when  specifying  boundary  conditions  for  numerical 
solutions  of  the  Navier-Stokes  equations  on  a  rectangular  domain. 

Our  objectives  are  to  formulate  and  solve  the  mathematical  system 
for  2-D,  unsteady,  standing  wave  disturbances  with  amplitudes  which 
vary  exponentially  in  the  streamwise  direction.  The  mean  flow  will  be 
represented  by  a  parallel-flow  boundary  layer.  The  first  evidence 
which  we  saw  of  the  existence  of  these  waves  was  in  an  analytical 
solution  of  spatially-decaying  vortices  in  the  freestream.  We  will 
summarize  those  results  in  the  next  section.  In  later  sections,  we 
will  extend  these  analytical  solutions  to  include  the  presence  of  a 
nearby  wall,  and  later  to  numerical  solutions  with  Falkner-Skan 
boundary  layers.  While  our  attention  here  is  focused  on  either 
oscillations  in  the  freestream  or  with  a  wall  (and  boundary  layer) 
present,  these  disturbances  should  arise  in  a  wide  variety  of  shear 
layers,  jets,  wakes,  etc. 
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1.2  STANDING  WAVES  AND  VORTICES  IN  A  UNIFORM  MEAN  FLOW 


The  Navier -Stokes  equations  can  be  separated  into  mean  and 
disturbance  quantities,  the  (time  or  ensemble)  average  taken,  those 
mean  equations  subtracted  from  those  equations  with  separated 
quantities,  and  simplified  for  a  parallel-flow  boundary  layer  with 
small-amplitude  disturbances.  With  appropriate  derivatives  in  the  x 
and  y  directions  taken,  a  fluctuating  vorticity  equation  can  be  found 
and  nondimensionalized.  For  a  uniform  mean  flow,  this  vorticity 
equation  is 


V 


t 


~e(f. 


XX 


V  - 


(1.3) 


where  ?  In  Refs.l,2,Rogler  and  Reshotko  sought  the  spatially 
decaying  solution  of  this  equation  subject  to  the  initial  condition  on 
vorticity 


f  ( *  2  *  stn(7rt/*yf)  s//i  vt 

and  the  far-downstream  condition  that  the  vorticity  vanishes 

f  O  as  x  -*■  co 


(1.4) 


(1.5) 


Conditions  on  the  velocity  field  are  imposed  later.  By  generalized 
Fourier  transforms,  a  solution  was  extracted  as 


x  +  irrt  mx-t7rt 


-  e 


] 


f  «  -<7rs<W(iry+y)  [e 

-  [/  -  (/^  */rzez  1 J// Ze 

or  equivalently 

?  =  -  2  77  5</?(  7T(/+yf )  e  ^,7r  7r£  (x-  c?) 


(1.6) 


(1.7) 


(l.B) 


This  result  is  a  slowly  decaying  wave  in  the  downstream  direction 
propagating  at  speed  c>  1. 

The  solution  of  interest  here  follows  from  Poisson's  equation 
which  links  the  vorticity  field  and  the  streamfunction 


=  f 

where  we  imposed  an  initial  condition  on  the  streamfunction 


(1.9) 
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(1.10) 

Again,  far  downstream,  the  streamfunction  vanishes 

</>  O  as  x  —  co 

(1.11) 

If  a  solution  of  form 

V  =  Jt*(?ry+S$)  €) 

(1.12) 

is  sought  of  Poisson’s  equation,  then  it  reduces  to 
&xx  -  7TZ0  =  Z7T  e  ^V* S(#  *(%*-?) 

(1.13) 

subject  to  the  conditions 

0(0\ and  &-~Oasx—<x> 

(1.14) 

The  complementary  solution  is 

ec(x ;  t,  e)  =  ct( t)  +  cg(t)  e"  W* 

and  it  is  this  solution  which  serves  as  the  initial  seed 

(1.15) 

for  this 

research.  We  shall  return  to  this  solution  after  we  have  written  the 
particular  integral,  obtained  by  the  method  of  variation  of  parameters. 
Combining  the  complementary  and  particular  solutions  and  applying  tne 
boundary  conditions,  then  the  disturbance  streamf unction  is 

t  =  £  s/'/!(vi/*v)[7rcz(t,  e)  e“ 

Term  (i )  Term  (Z) 

where  the  following  are  defined  ^ 

7TCz{t^)  =  \rfzC0S7tt  -f*  ] 

♦  [  4  cos  f)  ntWf)Z*£  f'- s<n  nt 

+  [('"  0  -]£  cos  77(/zX- ?)][(/- f)Z+  //]  Cl. 17a, b) 

Term  (1)  represents  a  near-grid  effect  which  is  an  irrotational, 
standing  wave  which  decays  exponentially  downstream.  The  second  term 
is  a  traveling  wave  and  represents  the  flowfield  of  the  convecting  and 
diffusing  vorticity  field.  The  streamlines  for  these  two  flowtields 
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are  plotted  in  Fig.  1.1 

The  standing  wave  was  "initiated"  by  the  initial  condition  imposed 
on  the  streamf unction.  It  would  have  been  possible  to  describe  initial 
conditions  in  which  that  standing  wave  vanished,  and  indeed  for  the 
case  studied,  this  wave  vanished  in  the  inviscid  limit.  With  other 
initial  conditions,  that  wave  would  remain  in  the  inviscid  limit. 


1.3  STANDING  WAVE  ARISING  IN  A  QUARTER-PLANE  PROBLEM 

In  Ref.  6,  a  more  complex  standing  wave  showed  up  in  an  inviscid, 
quarter-plane  problem  of  an  array  of  vortices  input  along  the  y-axis 
with  a  plate  positioned  along  the  x-axis  as  shown  in  Figure  1.2.  The 
mathematical  system  solved  was 


£  *  O  =4  f  =  -Zir  5C/?7T(x-t)  5 in (ny •<-(/')  (1.1B) 

vH  =  f  (1.19) 

f  -  ^  (1.20, 

?(a>  =  £  f'rtfry+y,)  sinzr(x-t)  (121) 

VZ^<<}=  O  (1.22) 


¥(<}=  sin7r(x-6)-  C  fory=otx>o 

¥(<)~0  /os'X=0i&>O 

<f{e)-+Oas  9  f/Jis  fancied  asx-*ee 

and  the  solution  was  extracted  by  a  half-range  Fourier  transform  in  the 
y-direction.  The  solution  for  the  streamf  unction  ¥OJ  was 

/  =  stnqc  ]  senir(x-t) 

Term  (At)  Term  (Az  ) 

+  f(z)  lirtty  sinvt 

Term  (s) 
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Figure  1.2  Disturbance  streamlines  at  different  times  as  vortices 
are  introduced  along  the  upstream  boundary  and  propagate 
downstream.  This  flowfield  consists  of  both  irrotational 
and  rotational  parts.  The  irrotational  part  is  composed 
of  both  standing  and  traveling  waves. 
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2  = 

F(i)  -  -  /may  |  C((?rz)cos  7ri  +  si(7n)  s//?7rz 

and  citarz)  and  sHtfz)  are  the  cosine  and  sine  integrals.  Term  (Al) 
denotes  the  flow  induced  by  the  vorticity  in  an  unbounded  flow,  with 
term  (A2)  representing  an  ir rotational,  traveling  wave  superimposed  and 
associated  with  the  impermeability  condition.  The  term  of  importance 
here  is  term  (B)  which  is  a  standing  wave,  and  oscillates  only  in  time. 
Far  downstream  of  the  grid,  it  vanishes. 


(1.25a,b> 


1.4  DECAYING  WAVES  WITH  SEMI-INFINITE  PLATES 

Decaying,  oscillating  waves  have  been  found  in  several  problems 
where  freestream  disturbances  interact  with  the  leading  edge  of  a 
semi-infinite  plate.  The  author  has  investigated  these  waves  for  many 
forms  of  freestream  disturbances,  including 

(1)  Arrays  of  square  vortices  (Refs.  1,7, and  6)  arrays  of  rectangular 
vortices  (Refs. 9  and  10),  and  oblique  plane  waves  of  vorticity  (Refs.  9 
and  10) . 

(2)  Vortex  sheets  (Refs.  9  and  10)  convecting  downstream  at  speed 
c  =  1.  This  is  interesting  case  which  clearly  illustrates  the  blockage 
caused  by  the  plate  and  differences  between  the  flows  above  and  below 
the  plate. 

(3)  Vortex  streets  (Refs,  9  and  10)  composed  of  potential  vortices 
in  a  stable  (Kantian)  configuration.  These  vortices  propagate  at  a 
speed  dictated  by  the  strength,  sense,  and  spacing  of  the  vortices. 
This  speed  can  be  less  than  or  greater  than  unity  (depending  on  the 
direction  of  rotation) ,  and  reduces  to  c  =  1  in  the  limit  of  vanishing 
circulation  or  infinite  spacing. 

(4)  Fluctuations  wholly  irrotational  (Refs.  9  and  10)  which  propagate 
at  arbitrary  speeds.  Many  examples  are  possible,  depending  on  the 
relative  strengths  and  phases  of  the  irrotational  waves.  While  this 
model  blows  up  exponentially  in  the  direction  normal  to  the  plate,  the 
model  is  useful  in  analyses  (and  perhaps  computations)  for  some  studies 
with  freestream  fluctuations. 

Despite  the  different  characters  and  speeds  of  these  freestream 
disturbances,  the  etreamf unctions  for  the  case  of  freestream 
disturbances  with  a  semi-infinite  plate  can  all  be  written  in  the 
unified  form 
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=  ?ra)(wf)  -  r(%yS)  (1.26) 

where  $^°4s  the  streamf unction,  for  the  freestream  disturbance  (if  the 
plate  were  not  present)  and  P(e*iB  the  streamfunction  of  an  irrotational 
flow  which  causes  the  impermeability  condition  to  be  satisfied.  We 
illustrate  this  second  streamfunction  in  Figure  1.3.  Interestingly, 
this  streamfunction  is  also  the  linearized  solution  of  a  semi-infinite, 
2-D  (i.e.,  ribbon- shaped)  snake  swimming  through  the  water. 


As  suggested  in  Figure  1.3  and  supported  by  the  theory  of  Refs.  9 
and  10,  the  streamfunction  has  a  simple  asymptotic  limit 
far-downstream  of  the  leading  edge 


a & =  <L(<> 

T  c  c  *  'asy'77/’  ■ 


(1.27) 


which  is  a  traveling,  irrotational  wave  which  decays  away  from  the 
plate.  If  (for  simplicity)  we  consider  the  quarter-plane  on  tne  top  of 
the  plate  beginning  at  the  leading  edge  and  extending  downstream,  then 
one  could  subtract  out  the  asymptotic  solution  from  The 
result  is  an  irrotational  flow  decaying  downstream  of  tne  leading  edge. 


,,  <t)  y/  (<)  _ 

*  ~  °  aSX  °°  <1.28> 

If  this  flow  (  ftr>-  )  were  Fourier-analyzed  in  tne 

y-direction,  then  one  would  obtain  the  spectrum  of  waves  which 
oscillate  neutrally  in  the  y-direction  and  decay  exponentially 
downstream.  These  spectral  components  are  related  to  tne  spectral 
waves  investigated  in  this  report.  Thus  we  believe  that  such  waves 
arise  for  a  wide  class  of  practical  freestream  disturbances  which 
intersect  with  the  leading  edges  of  plates. 


1.5  DECAYING  AND  GROWING  WAVES  WITH  FINITE-LENGTH  PLATES 

If  the  plate  has  finite  length,  with  both  leading  and  trailing 
edges,  then  the  resultant  flow  has  some  features  similar  to  those 
summarized  in  the  previous  section,  although  made  more  complicated  by 
the  trailing  edge.  The  flow  also  has  some  new  features  associated  with 
the  oscillating  circulation  about  the  plate  and  the  vortex  sheet 
emitted  from  the  trailing  edge  and  convecting  downstream. 

For  the  purpose  of  identifying  physically-realizable  situations 
where  standing  waves  grow  or  decay  in  the  streamwise  direction,  we 
offer  the  following  observations  based  on  the  analyses  of  Ref.  7: 

(1)  The  irrotational  flow  associated  with  the  bound  circulation 
oscillates  in  time  and  has  streamlines  as  plotted  in  Fig.  1.4.  This 
flow  is  irrotational,  oscillates  sinusoidally  in  time,  decays  from  the 
leading  edge  to  the  midspan,  and  grows  from  the  midspan  to  the  trailing 
edge. 
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Figure  1.3  The  streamlines  for  an  irrotational  flow  which  arises 
when  a  semi-infinite  plate  encounters  a  traveling-wave 
freestream  disturbance.  The  leading  edge  is  at  the  origin. 
Downstream  of  the  leading  edge,  this  flowtield  consists 
of  a  superposition  of  standing  waves  and  a  traveling  wave. 
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Figure  1.4  Streamlines  for  the  oscillating,  bound  circulation  about 
a  finite-length  flat  plate  (Ref .7).  This  flow  decelerates 
from  the  leading  edge  to  the  mid-span,  and  accelerates 
from  the  mid-span  to  the  trailing  edge.  This  flowtield  is 
irrotational  and  oscillates  sinusoidally  in  time.  It 
is  one  of  the  fluctuations  arising  when  an  airfoil 
encounters  a  traveling-wave  freestream  disturbance. 
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(2)  The  trailing  vortex  sheet  induces  an  oscillating,  irrotational  flow 
about  the  plate  as  shown  in  Fig.  1.5. 

(3)  To  satisfy  impermeability,  another  irrotational  wave  is  produced 
which  is  related  to  the  wave  discussed  in  the  previous  section.  If  tne 
plate  is  very  long,  then  this  wave  has  the  character  of  standing  waves 
near  the  leading  and  trailing  edges.  A  traveling  wave  is  also  present 
everywhere  along  the  plate.  The  streamline  pattern  associated  with 
this  flow  has  not  been  drawn;  the  explicit  solution  is  presented  in 
Ref.  7. 

In  summary  for  freestream  disturbances  interacting  with 
finite-length  plates,  growing  and  decaying  standing  waves  are  present 
simultaneously.  Several  mechanisms  exist  to  produce  standing  waves. 
While  a  unified  theory  is  available,  no  attempt  has  been  made  to 
compare  the  relative  strengths  of  these  various  waves. 


1.6  RELATED  STUDIES  OF  FOURIER-LAPLACE  AND  POTENTIAL  FLUCTUATIONS 

This  study  is  concerned  with  the  solutions  of  the  incompressible 
parallel-flow  equations  for  a  boundary  layer  on  an  impermeable  plate. 
To  help  provide  some  perspective  for  the  various  Fourier-Laplace 
solutions,  the  following  references  are  cited.  Of  course,  many  other 
investigators  have  made  significant  contributions  to  analyze, 
calculate,  and  measure  these  oscillations,  and  to  extend  and  apply  them 
to  many  other  problems. 

Tollmien  (Ref.  11)  first  obtained  neutral  eigenvalues  and  a 
critical  Reynolds  number  for  the  fundamental  stability  wave.  Jordinson 
(Ref.  36),  Mack  (Ref.  12),  Antar  and  Benek  (Ref.  37),  Corner,  Houston, 
and  Ross  (Ref.  13)  obtained  the  (temporal  and/or  spatial)  eigenvalues 
for  the  higher  discrete  modes  with  a  Blasius  layer.  Rogler  and 
Reshotko  (Refs.  1,2,3),  the  present  author  (Refs. 15  and  39),  Murdock 
(Ref.  40),  Salwen  and  Grosch  (Refs.  41,14),  and  Ellinwood  (Ref.  42) 
obtained  solutions  with  downstream  propagating  vortical  fluctuations  in 
the  freestream.  The  forced  response  of  a  boundary  layer  to  a  Karman 
vortex  sheet  was  obtained  by  Rogler  (Ref.  38).  Rogler  and  Reshotko 
(Ref.  2),  Rogler  and  Tsuge  (Refs.  17),  and  the  present  author  (Ref.  27) 
obtained  solutions  for  waves  travelling  in  the  upstream  direction. 
Those  authors  (Ref. 16)  also  obtained  solutions  for  the  standing  waves, 
with  details  summarized  in  the  present  report.  In  Ref.  17,  spatial 
solutions  for  longitudinal  vortices  in  a  boundary  layer  were 
summarized.  Related  temporal  solutions  were  obtained  by  Rogler  and 
Reshotko  (Ref.  18)  and  Hultgren  and  Gustavsson  (Ref.  19).  The 
boundary- value  and  initial-value  problems  will  be  reviewed  in  a  later 
report. 

For  temporal  solutions  of  the  Orr-Sommerfeld  equation  in  a 
channel,  Mack  (Ref.  12)  presented  numerical  evidence  of  how  the 
discrete  solutions  for  a  channel  carry  over  to  the  solutions  for  a 
semi-infinite  domain  ( 0^y<“0 .  There  should  be  a  link  between  the 
standing  waves  of  the  present  report  (ase-^0)  and  oscillations  in  a 
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5  Streamlines  of  the  . flow  induced  by  the  vortex  sheet 
emitted  from  the  trailing  edge  of  a  finite-length  plate 
(Ref.  7)  .  Only  the  top  half  has  been  drawn.  Except  at  tne 
vortex  sheet  itself,  the  flow  is  irrotational ,  oscillates 
sinusoidally  in  time,  and  has  characteristics  of  both 
traveling  and  standing  waves. 


»-EB-til003V 


AEDC-TR-83-4 


channel  (as  the  channel  width  —+■»  and  the  x-wavenumber  -*>  0) .  The 
neutral  oscillation  v=0,  u=f  (y)exp(-i<»t)  should  be  a  limiting  solution 
for  both  cases. 

The  importance  of  these  Fourier-Laplace  solutions  are  that 

(1)  they  can  be  superimposed  to  represent  more  general  disturbances, 

(2)  the  separate  solutions  contain  many  of  the  features  arising  through 
instabilities,  walls,  viscosity,  details  of  the  mean  shearing  layer, 
roughness,  freestream  disturbances,  etc.,  and  (3)  the  Fourier-Laplace 
tools  lead  to  a  powerful  and  flexible  framework  for  the  study  of 
transition  and  mixing.  This  framework  is  important  because  progress  is 
often  based  upon  generations  of  work  by  many  investigators,  a  diversity 
of  (analytical,  numerical,  and  experimental)  tools,  and  engineering  and 
scientific  problems  of  many  types.  The  usefulness  of  stability  theory 
has  been  firmly  established.  Now,  the  potential  exists  for  a  broader 
theory,  which  includes  the  normal  modes,  but  also  incorporates  some 
features  of  freestream  disturbances,  surface  waviness  and  roughness, 
surface  vibration,  and  forcing. 

The  reader  is  reminded,  however,  that  not  all  monochromatic 
Fourier-Laplace  disturbances  in  the  freestream  yield  a  monochromatic 
response  in  a  parallel-flow  boundary  layer,  nor  even  in  a  uniform  mean 
flow  near  a  wall.  In  the  temporal  initial-value  problems  which  have 
been  worked,  the  wavenumber  of  the  freestream  disturbance  is  preserved, 
but  the  frequency  is  transformed  into  a  spectrum.  In  the  spatial 
boundary-value  problems  carried  out  by  the  author  for  an  inviscid  flow, 
the  frequency  of  the  freestream  disturbance  is  preserved,  but  an 
integral  spectrum  and/or  a  series  of  wavenumbers  result. 

The  work  by  Bechert  and  Michel  (Ref.  20)  is  an  example  of  a 
potential  fluctuation  in  the  form  of  an  oscillating  source/sink  used  to 
excite  a  step-function  free-shear  layer  in  an  inviscid  analysis.  The 
experimental  investigation  of  Dovgal,  Kozlov,  and  Levchenko  (Ref.  21) 
is  believed  to  include,  amongst  the  various  fluctuations,  the  standing 
wave.  In  other  experiments,  the  author  cannot  distinguish  whether 
compressibility  is  an  essential  feature  of  the  input  disturbance,  or  if 
the  input  disturbance  is  partly  composed  of  incompressible  standing 
waves.  Other  incompressible  feedback  problems  with  free-shear  layers 
emitted  from  trailing  edges,  and  the  incompressible  coupling  between 
the  trailing  edge  and  a  leading  edge  positioned  further  downstream  also 
have  some  elements  of  the  standing  waves  which  are  investigated  here. 
Additionally,  some  of  the  incompressible,  potential  fluctuations  in 
unsteady  airfoil  theory  also  have  certain  characteristics  of  the 
standing  waves,  as  discussed  in  Section  1.5. 

Hence  the  present  work  might  be  viewed  as  a  contribution  toward 
bringing  these  additional  features  into  an  Orr-Sommerfeld  description 
of  unsteadiness  in  a  shearing  layer.  As  discussed  by  Reshotko 
(Ref.  22),  a  mathematical  and  physical  understanding  of  tne  processes 
at  the  beginning  of  transition  "would  be  most  welcome." 

In  the  next  chapter,  the  mathematical  system  describing  the 
standing  waves  in  viscous,  parallel-flow  boundary  layers  is  formulated. 
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2.0  FORMULATION  OF  THE  MATHEMATICAL  SYSTEM 
2.1  THE  ORR-SOMMERFELD  EQUATION  AND  BOUNDARY  CONDITIONS 

The  linearized,  parallel-flow,  viscous  momentum  equations  can  be 
combined  into  the  following  fourth-order  partial  differential  equation 
in  terms  of  the  fluctuating  velocity,  v(x,y,t) ,  in  the  y-direction  (and 
normal  to  a  plate) 


The  inverse  Reynolds  number  based  upon  the  freestream  velocity  and 
characteristic  thickness  of  the  boundary  layer  is  e  =  *>71^5.  For 
Falkner-Skan  boundary  layers,  the  characteristic  thickness  is  taken  as 


where  *  frs/(  z~  fi *s)  (2.2a,b> 

is  the  vertex  angle  of  the  wedge  along  which  the  boundary  layer 
develops.  The  lengths  in  eqn.  2.1  have  been  nondimensionalized  against 
&  and  the  time  has  been  nondimensionalized  against  S/EJ^.  We  now  seek 
solutions  of  form 


v(*>y>  *)  ~  9(y)e 


MX  -  /cot 


(2.3) 


where  the  x-wavenumber  (times  i)  for  the  cases  of  interest  here  are 


=  7?  or  '*=-<> p  (2.4) 

where  8  is  pure  real,  and  where  the  coefficient  r  has  values  +1  and  -1 
corresponding  to  the  cases 


r  =  +1  for  a  wave  growing  as  exp(+&t) 

r  =  -1  for  a  wave  decaying  as  exp(-px)  (2.5a,b> 

The  x-wavenumber  is  pure  imaginary.  Substituting  eqn.  2.3  into  2.1, 
then  after  some  rearrangement,  the  Orr-Sommerfeld  equation  results 


\{v~  -  q,  t -o 

or  equivalently 

\(o-  -fa  (*+?)'} f = ° 

The  phase  speed  is  pure  imaginary  for  the  standing  wave  solutions 
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C  H  a>/<*  = 


(2.7) 


Either  equation  is  subject  to  impermeability  and  no-slip  at  the  wail 

=  Pf  =  O  af  y-O  (2.8) 

When  U  =  1  outside  the  boundary  layer ,  then  the  system  reduces  to 


=  o 


(2.9) 

r  °  j 

±  l&U 

with  solutions  e  ' ^  or  cosay  and  singy  (which  oscillate  neutrally) 
±mu  ' 

and  €  (which  oscillate  and  grow  or  decay).  (2.10arb) 


where  m  is  the  root  of  the  quadratic  equation  with  positive  real  part 

=  O 

(2.11) 


<  p  '  1  r  ' 

,  I fz 

or  /7?  -  t  with  r  =  +1  for  exp(+px)  waves  (2.12) 

The  exterior  solution  is 

(e>(?)  =  de  y+3e  Ccos^y  +  (2.13) 

where  B  =  0  in  the  present  problem  and  we  normalize  by  setting  C  =  1. 
This  normalization  is  not  always  appropriate;  see  Appendix  A.l  for  an 
exception.  We  require  that  the  exterior  solution  and  the  first  three 
derivatives  agree  with  the  interior  solution  and  the  corresponding 
derivatives  at  the  "edge"  of  the  boundary  layer,  y=y 


/o~  *7*  Ot  /,  2}  3 


(2.14) 


Our  objective  now  is  to  solve  the  mathematical  system  composed  of 
the  Orr-Sommerfeld  equation  (eqn.  2.6a,b),  the  two  wall  conditions 
(eqn.  2.8),  and  the  four  conditions  at  the  boundary  layer  edge 
(eqns.  2.13  and  2.14).  The  exponents  m  and  the  other  parameters  for 
the  cases  of  decaying  and  growing  waves  are  given  in  Table  1  as  a 
typical  example.  Note  that  the  exponent  m  for  the  decaying  case  is 
nearly  i  times  m  for  the  growing  case.  The  slight  difference  arises 
from  the  term  -sz  in  the  square  root. 
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2.2  COMPLEX  AMPLITUDES  FOR  LONGITUDINAL  VELOCITY,  VORTICITY, 
STREAMF UNCTION,  AND  PRESSURE 


Let  the  velocity  fluctuations  in  the  x  and  y  directions  be 
represented  as 

i'OhrwV**-'* 

Substituting  into  the  continuity  equation,  ux+  ^,=0,  then  (<*/  +  fy-=  O 
and  the  amplitude  of  the  longitudinal  velocity  fluctuation  is  related 
to  the  amplitude  of  the  normal  velocity  fluctuation  by 


S(yJ  =  <<%(?)/« 


(2.16) 

The  fluctuating  vorticity  is  f~  Vx-Wy  •  Letting  f  , 
then  the  vorticity  amplitude  Is  related  to  the  two  velocity  amplitudes 
by 


(2.17) 


Letting  ^  ^ 


since  the  normal  velocity  and 
streamfunction  are'  related  by  v  ■*  ,  then  the  amplitude  of  tne 

st reamf unction  is 


¥(</)  = 


(2.18) 


The  pressure  disturbance  can  be  found  from  the  linearized, 
parallel-flow,  x-momentum  equation 


tSc/V  —  °  +  *  ^yy) 


(2.19) 


Introducing  the  dimensionless  variables  u-*qu,  v—qv,  t-*£t/U  ,  U-^l^U, 
p -*oqUp  ,  x-*<£x,  and  y-**Sy,  where  q  is  the  characteristic  disturbance 

velocity,  then 


II 

-£■ 

*&(*****») 

(2.20) 

Letting 

P  * 

TT(y)e 

(2.21) 

then 

II 

1 

Yf 

(2.22) 

Thus  knowing  the  amplitudes  of  the  velocity  fluctuations  and  the 

mean 

velocity 

profile. 

then 

the  pressure  fluctuation  can  be  found. 

In  a 

manner  analogous  to 

the 

streamlines,  the  pressure  isobars  can  be 

calculated  and  plotted. 
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2.3  ROOT-MEAN-SQUARE  QUANTITIES  AND  OTHER  CORRELATIONS 


Let  a  disturbance  quantity be  represented  as 

G(°(xt^)  =  iff 

=  6>r0){y)cos( *X-£ot)  -  (<*X-  (tit )  (2.23) 

Similarly#  a  second  disturbance  quantity  ^^can  be  represented  as 

e™(y)COs(°<X-wt)~  &<*fc)s/s?frx-&?)  (2.24) 
The  product  of  these  two  (real)  quantities  is 

o>  (a)  -  (/)  (Z)  .2 

G  G  -  <2  COS  (*X-6U?)  +  &■  5<v?  (- Oi-£& ?) 

-(G*'W*+  ef'ty**)  ros(*X-wt)s<n(«X-aft)  (2.25) 

By  the  mean  of  a  variable  in  a  spatial  analysis#  we  refer  to  the  time 
average  over  one  time  period  T  =  2Tf/co  of  that  variable 

T 


&(°  e(z) 


(2.26) 


Carrying  out  the  integrals  for  & f^*;the  contributions  arise  only  from 
the  squared  quantities,  and  the  result  is 


(/)  (z)  /  /  _  (O  <*)  6*7) 

*  *  -  4  (°r  8r  *  6i  *<  J 


(2.27) 


The  autocorrelations  are  thus 


&)¥* - u‘~  ?)M/A 

(^)/z  =  O/zf 
(7~)/z  =  [(*•♦ 

The  Reynolds  stress  is  —  t/V  -  -  (£  ^  fyj/z 

The  energy  production  is  —  6/ V 


(2.28) 


(2.29) 


(2.30) 
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The  {instantaneous)  kinetic  energy  per  unit  mass  is 

■j(#Z+V*)  =  J  (o*+ Z Da  *  v* jrf- +  v3)  (2  31) 

The  mean  of  this  quantity  is 

j.(oz  +  Z7*  7r) 

2  (  /  (2.32) 

which  has  a  contribution  from  the  mean  flow  and _ another  contribution 

from  the  unsteady  flow.  The  quantity  (u*+va)/2  is  plotted  in  the 
following  chapter,  along  with  the  other  variables  introduced  above. 


Table  1 

Example  Numerical  Value  of  the  Exponent  m 


exp(-(3x) 

mode 

P 

0.5 

cu 

0.5 

1000 

* 

0.5i 

/7i 

10.1743201  -i24. 5716664 

exp(+jjx) 

mode 

0.5 

0.5 

1000 

-0.5i 

24.5625)806  -ilO. 1779179 
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3.0  NUMERICAL  SOLUTIONS  FOR  STANDING  WAVES 


3.1  SUMMARY  OF  THE  NUMERICAL  TECHNIQUE 

To  obtain  the  coefficients  U(y)  and  Uyy(y)  of  the  mean  flow, 
numerical  solutions  of  the  Falkner-Skan  equation  are  obtained  by  the 
method  of  Nachtsheim  and  Swigert  (Ref*  23)  as  modified  and  summarized 
in  Ref.  3.  The  equation  is  solved  by  shooting  from  the  wall  outward, 
using  4th-order  Runge-Kutta  integration  in  double-precision  arithmetic. 
The  exterior  boundary  condition,  R>-*'1  as  ,  is  not  specified  at 

the  boundary  layer  edge,  y  *=  ,  but  ^rather  the  error  is  minimized  in  a 

least  squares  sense. 

The  Orr-Sommerfeld  system  of  equations  (2.6b,  2.8,  2.13,  2.14)  is 
solved  by  an  expansion  in  Chebyshev  polynomials 

<P(V  =YL°mT,„(?) 

(3.1) 

with  polynomials  defined  as 

cos[m(arccos?y\  (3.2) 

and  where  the  variable  ^  is  related  to  the  dimensionless  y  by  the 
linear  transformation 


A  set  of  linear,  algebraic  equations  results  from  the  Orr-Sommerfeld 
equation,  the  two  wall  boundary  conditions,  and  tne  four  matching 
conditions  at  the  boundary  layer  edge.  The  matrix  of  coefficients 
related  to  this  system  is  reduced  by  Gauss-Jordan  elimination  with 
maximum  pivoting  by  column  and  rows.  With  the  coefficients  known,  the 
expansion  3.1  can  be  carried  out  for  ^(y),  and  identities  used  to 
obtain  the  y— derivatives  of  ^(y) •  With  and  its  derivatives 
known,  the  longitudinal  velocity,  vorticity,  and  vorticity  production 
are  calculated.  From  the  x-momentum  equation,  the  pressure 
fluctuations  are  found,  and  then  the  correlations  are  calculated  for 
rras  quantities,  including  Reynolds  stress,  kinetic  energy,  etc.  The 
streamlines,  iso-vorticity  contours,  pressure  isobars,  and  equi-value 
contours  of  the  vorticity  production  term  are  drawn. 

The  numerical  checks  included  numerical  experiments  to  determine 
the  effects  of  number  of  terms  in  the  expansion,  precision  of 
arithmetic  for  the  matrix  reduction,  and  influence  of  the  y-value  where 
the  edge  boundary  conditions  were  imposed.  Comparisons  were  made 
between  analytical  and  numerical  solutions  for  the  cases  of  (l)a 
uniform  mean  flow,  and  (2)  the  viscous  sublayer.  The  same  computer 
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program  used  to  obtain  solutions  for  the  standing  waves ,  except  witn 
different  patching  conditions  at  the  boundary  layer  edge;  was  also  used 
to  recover  the  Tollmien-Schlicting  waves. 

To  assist  the  reader  in  comparing  the  figures  with  different 
values  of  the  parameters ,  different  wall  boundary  conditions,  and  tne 
two  cases  of  growing  or  decaying  waves.  Table  2  has  been  prepared. 


3.2  NUMERICAL  RESULTS  FOR  DECAYING  WAVES  IN  FALKNER-SKAN  BOUNDARY 
LAYERS 


3.2.1  Blasius  Case  with  Decaying  Waves 

The  velocity  profile  and  derivatives  for  the  Blasius  boundary 
layer  are  plotted  in  Figure  3.1  with  tne  same  y-scale  as  tne  following 
figures  for  amplitudes,  streamlines,  etc,,  except  for  the  sublayer 
plots  where  the  y-scale  is  stretched. 

To  aid  in  the  interpretation  of  these  calculations,  since 

—  fix—  -fix 

V  -  (<Pr  +  e'4>{  )  e  “  =  [  <Pr(y)<-osa)t  +  J  e  “ 

(3.4) 

then  for  cut  -  O  ,  v(x>y,0)-  ^(y)? 

6ot^?r/2y  V(xtytt)  =  <pt(y)e~?x 
a>t  =  7r,  v(x,y,  ?)  = 

^  =  3^/2,  V{x,yft)  =  -<$-(y)e~P*  (3.5a-d) 

Hence  the  real  and  imaginary  parts  of  the  complex  amplitude  are  the 
velocity  profiles  at  the  times  oot=Q  and  7?/2  respectively.  With  a 
reversal  of  sign,  they  are  the  velocity  protiles  at  the 
times  wt.=  7T  and  3#V2.  At  <^t=2 7f ,  the  flowfield  at  time  t=0  is 
repeated  in  this  sinusoidal  oscillation.  Analogous  interpretations 
also  apply  for  the  other  complex  amplitudes. 

Figures  3.2  through  3.19  are  plots  for  a  decaying  wave  with 
y-wavenumber  and  frequency,  p  =<^=  0.5,  and  Reynolds  number  Rj=1000.  In 
Figure  3.2  for  the  normal  velocity,  note  that  the  impermeability 
condition  is  satisfied  at  the  wall  and  note  the  oscillatory  behavior  of 
the  solution.  Far  from  the  boundary 

<?(</)  —  StKfy  +  p  as ao  (3>6) 

Also  note,  based  on  the  above  interpretation  of  the  real  and  imaginary 
parts,  that  the  flow  structure  shifts  up  as  the  time  increases  from  0 
to  W/7.<u,  The  structure  nearest  to  the  wall  expands  outward.  To  see 
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Table  : 

2 .  Tabular-List 

of  Figures 

for  Chapter  3 

-ax  +bx 

e  '  or  el  wave 

e”Px 

e-px 

e-P* 

Wall  BC 

no-slip 

- >■ 

no-shear 
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0.5i 

0.5i 

0.5i 

-0.5i 

0.5i 

frequency, co 
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0.5 
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Figures  illustrating  the  effects  of  frequency,  Reynolds  number,  and  y-wavenumber 
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this,  observe  that  the  node  where  <fy=0  shifts  upward  to  the  node 
where  ^=0. 

In  Figure  3.3  for  the  longitudinal  velocity,  note  that  the  no-slip 
condition  is  satisfied  and  again  note  the  general  oscillatory  behavior 
further  away.  The  fluctuating  velocity  inside  tne  boundary  layer  is 
much  larger  than  the  velocity  in  the  boundary  layer  with  decaying, 
vortical  structures  in  the  freestream. 

In  Figure  3.3,  also  note  the  thin  viscous  sublayer  near  the  wall, 
with  a  velocity  overshoot  at  the  sublayer  edge.  In  Figure  3.4,  the 
ordinate  is  stretched  and  the  longitudinal  velocity  is  plotted  for  the 
sublayer  region.  The  analytical  solution  for  the  sublayer,  as 
presented  in  Appendix  B,  is  also  plotted  in  Figure  3.4.  This  solution 
is 


where 


(3.7) 


and  5J®  -  complex  constant  representing  the  value  of  f  at  the 
sublayer  "edge". 

The  streamlines  for  the  fluctuating  flow  at  time  zero  are  plotted 
in  Figure  3.5.  Note  the  oscillatory  pattern  in  the  y-direction  and  tne 
effects  of  the  exponential  decay  in  the  streamwise  direction, 
exp(-jSx).  A  16mm  movie  of  the  standing  wave  has  been  prepared  using  a 
sequence  of  64  streamline  patterns  representing  the  flowfield  at  64 
times  during  a  half -cycle. 

The  fluctuating  vorticity,  f  =vx-Uy,  is  plotted  in  Figure  3.6, 
which  shows  a  large  spike  of  vorticity  near  the  wall  associated  with 
the  viscous  sublayer.  This  vorticity  mainly  arises  from  the 
derivative  -u*,  while  the  derivative  v*  is  zero  or  small.  The  large 
magnitude  of  the  derivative  a„  in  the  sublayer  can  be  seen  in  Figures 
3.3  and  3.4. 


The  abscissa  in  Figure  3.6  was  stretched  in  Figure  3.7  to  better 
display  the  layer  of  vorticity  appearing  about  midway  in  the  boundary 
layer;  the  much  larger  spike  of  vorticity  in  the  sublayer  is  not 
completely  plotted  in  Figure  3.7.  Note  in  either  Figure  3.6  or  3.7 
that  far  above  the  boundary  layer  into  the  freestream,  the  flow  is 
irrotational.  The  iso-vorticity  contours  are  plotted  in  Figure  3.8 
which  shows  the  highly  skewed  nature  of  the  vorticity  waves  in  the 
sublayer.  The  single  contour  toward  the  center  of  the  boundary  layer 
represents  the  bulge  in  vorticity  near  y  =  2  in  Figure  3.7. 

The  fluctuating  vorticity  positioned  approximately  midway  in  the 
boundary  layer  has  interesting  properties  which  we  will  describe  by 
posing  two  mathematical  systems  for  inviscid  flows; 
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Figure  3.3  The  variation  of  the  longitudinal  velocity  fluctuation 
with  distance  from  the  wall  for  a  decaying  wave# 
exp(-jflx)  ,  for  /3  =<&=  0.5  and  Rj-1000  with  a  Blasius 
boundary  layer. 
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Figure  3.4  The  longitudinal  velocity  fluctuation  in  the  sublayer 
for  a  decaying  wave,  exp(-ySx)  ,  for  e  =eO=  0.5  and  R5=1000 
with  a  Blasius  boundary  layer. 
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Figure  3.5  The  streamlines  of  the  unsteady  flowfield  for  a  decaying 
wave,  exp(-^x)  ,  for  0.5  and  R^=1000  with  a  Blasius 

boundary  layer. 
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Figure  3.6  The  variation  of  the  vorticity  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave,  expf-^x) ,  for 
p=u>=  0.5  and  R^=1000  with  a  Blasius  boundary  layer.  This 
figure  shows  the  vorticity  in  the  sublayer.  Figure  3.7 
shows  the  vorticity  in  the  middle  of  the  boundary  layer. 
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Complex  amplitude  of  the  vorticity  fluctuation,  z(y) 


Figure  3.7  The  variation  of  the  vorticity  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave,  exp(-jsx)  ,  for 
/S  =  ry  =  0.5  and  Rj-1000  with  a  Blasius  boundary  layer. 
The  abscissa  has  been  stretched  to  show  the  details  of 
the  vorticity  in  the  center  of  the  boundary  layer. 
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Figure  3.8  The  iso-vorticity  contours  of  the  unsteady  flowfield  for 
a  decaying  wave,  exp(-^x) ,  for  /S=<w=  0.5  and  1^=1000 
with  a  Blasius  boundary  layer.  Fifty-seven  contours  were 
plotted  with  unequal  increments  of  vorticity. 
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(1)  If  there  were  no  vorticity  production  terra,  the  relevant 
vorticity  equation  would  be 


5  *  v(</)  fx  -  O 


(3.8) 


i.e.,  the  vorticity  at  any  station  x,y  at  time  t  is  the 
vorticity  at  x=0  at  the  earlier  time  t-x/U(y).  The 
vorticity  would  merely  convect  downstream  with  the  local  mean 
velocity,  U(y).  Since  the  vorticity  at  x=0  oscillates  in 
time 


(3.9) 


the  solution  further  downstream  would  be 


10) 


which  is  a  travelling  wave,  in  disagreement  with  the  assumed 
standing  wave  solution  of  form 

_fix-  UO<? 

T  =  Z(y)e  £3.11) 

One  may  wonder  then,  for  a  standing  wave,  what  causes  the 
vorticity  to  decay  exponentially  in  the  streamwise  direction 
(at  one  instant  of  time)?  To  an  observor  travelling  with  the 
mean  flow,  U(y),  how  can  the  vorticity  oscillate? 


(2)  We  believe  that-  one  mechanism  which  influences  the 
vorticity  to  behave  as  a  standing  wave  is  the  vorticity 
production  term,  -vC Ly ,  which  is  a  linearized  version  of  the 
convection  term  Vj£  in  the  vorticity  equation.  The 
linearized,  inviscid  equation  for  fluctuating  vorticity  is 

+  V%  =  °  (3.12) 

rate  of  change  of  vorticity  rate  of 
as  seen  by  an  observor  moving  production 
with  the  local  mean  flow 

For  high  Reynolds  number,  small  y-wavenumber  conditions,  we 
suspect  that  fluctuating  vorticity  is  created  in  the  central 
region  of  the  boundary  layer,  as  described  by  the  production 
term,  in  such  a  manner  that  a  standing  wave  is  possible.  The 
magnitude  of  the  vorticity  which  is  produced  depends  on  the 
time  available,  and  hence  the  frequency.  While  this  process 
occurs  in  the  central  region  of  the  boundary  layer,  diffusion 
is  possible  in  the  sublayer.  The  relative  roles  of  the 
production  term  and  the  two  diffusion  terms  (in  the  viscous 
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equation)  could  be  determined  by  plotting  the  terms 

separately. 

The  complex  amplitude  of  the  vorticity  production  is  plotted  in 
Figure  3.9.  The  production  vanishes  at  the  wall  and  in  the  freestream. 
The  contours  of  equi-value  of  the  production,-v(x,y,t) Uj^ty)  at  time 
zero  are  plotted  in  Figure  3.10. 

The  diffusion  term  in  the  streamwise  direction 

is  merely  6  &Z  times  the  vorticity  itself.  The  other  diffusion  term, 
€  ,  has  hot  been  plotted. 

The  smoothly  oscillating  amplitudes  57"  ( y)  =  7Tr+*7t re-  of  tne 
fluctuating  pressure  are  plotted  in  Figure  3.11.  Note  the  relative 
constancy  of  the  pressure  across  the  sublayer,  but  clearly  the  pressure 
varies  across  the  mean  boundary  layer  itself.  The  corresponding 
pressure  isobars  are  plotted  in  Figure  3.12. 

The  rms  quantities  for  the  velocities,  vorticity,  and  pressure  are 
plotted  in  Figures  3.13-3.16.  The  oscillatory  behavior  in  the 
direction  normal  to  the  wall  is  apparent  in  these  figures.  The  viscous 
sublayer  is  shown  clearly  by  the  longitudinal  velocity  in  Figure  3.14, 
and  by  the  spike  of  vorticity  near  the  wall  in  Figure  3.15. 

The  Reynolds  stress,  -uv,  and  the  kinetic  energy  production, 
-uvdU/dy ,  are  plotted  in  Figures  3.17  and  3.18.  The  Reynolds  stress 
oscillates  into  the  freestream.  The  energy  production  is  confined  to 
the  boundary  layer,  of  course,  where  it  is  mainly  negative  although 
having  a  small  positive  region.  It  vanishes  at  the  wall  where  v=0  and 
in  the  freestream  where  =0 . 

The  averaged  kinetic  energy  tuz+vz)/2  of  the  fluctuation  flow  is 
plotted  in  Figure  3.19  which  clearly  shows  the  sublayer  and  overshoot 
at  the  sublayer  edge.  In  the  freestream,  the  kinetic  energy  is  uniform 
in  the  y-direction;  it  decays  as  exp(-2^x)  in  the  streamwise  direction. 


3.2.2  Favorable  find  Adverse  Pressure  Gradient  Gases  with  Decaying  Waves 

Figures  3.20-3.26  are  results  for  the  decaying  wave  with  a 
favorable  pressure  gradient  boundary  layer.  Figures  3.27-3.33  are 
results  with  an  adverse  pressure  gradient  boundary  layer.  The 
corresponding  figure  numbers  for  these  two  cases  and  the  Blasius  case 
are  tabulated  below  for  reference. 
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Figure  3.9  The  variation  of  the  vorticity  production  with  distance 
from  the  wall  for  a  decaying  wave,  exp(-^x) ,  for 
p  0,5  and  R5=1000,  and  a  Blasius  boundary  layer. 
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Figure  3.10  The  equi-value  contours  of  the  production  of  fluctuating 
vorticity  for  a  decaying  wave,  exp(-^x)  ,  for  p  =  cu=  0.5, 
Rj=1000  at  time  t  =  0  with  a  Blasius  boundary  layer. 
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Figure  3.11  The  variation  of  the  pressure  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave,  exp(-ysx) ,  for 
a  o.5  and  Rg-1000  with  a  Blasius  boundary  layer. 
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Figure  3.13  The  variation  of  the  rms  normal  velocity  fluctuation  with 
distance  from  the  wall  for  a  decaying  wave,  expf-^x) ,  for 
0,5  and  R^=1000,  and  a  Blasius  boundary  layer. 
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Figure  3.14  The  variation  of  the  rms  longitudinal  velocity 
fluctuation  with  distance  from  the  wall  for  a  decaying 
wave,  exp(-ax) ,  for  0.5  and  1^=1000  with  a  Blasius 

boundary  layer. 
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Figure  3.15  The  variation  of  the  rms  vorticity  fluctuation  with 
distance  from  the  wall  for  a  decaying  wave,  exp(-/3x)  ,  for 
p  =cv=  0.5  and  R&»1000  with  a  Blasius  boundary  layer. 
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Figure  3.17  The  variation  of  the  Reynolds  stress  with  distance 
the  wall  for  a  decaying  wave  for  p  =<2c/=  0.5  and 
with  a  Blasius  boundary  layer. 
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Figure  3.19  The  variation  of  the  kinetic  energy  with  distance  from 
the  wall  for  a  decaying  wave,  expC-^x)  ,  for  0.5 

and  R&=1000  with  a  Blasius  boundary  layer. 
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A  plot-by-plot  comparison  o£  the  three  cases  shows  that  the  solutions 
for  the  standing  wave  depend  on  the  Falkner-Skan  parameter.  The 
maximum  values  and  the  y-positions  for  those  maximum  values  depend  on 
the  mean  profile. 

The  general  trend  is  that  the  curves  are  shifted  upward  aa  the 
Falkner-Skan  parameter  decreases.  The  energy  production  acquires  a 
positive  region  near  the  boundary  layer  edge  as  the  Falkner-Skan 
parameter  takes  on  negative  values. 

Based  upon  the  idea  that  one  effect  of  surface  roughness  is  to 
alter  the  mean  velocity  profile,  then  our  calculations  witn  different 
Falkner-Skan  parameters  indicates  that  the  surface  roughness  would 
influence  the  standing  waves. 


3.2.3  Effects  qL  Frequency 

As  shown  in  Figures  3.34  and  3.35,  the  frequency  significantly 
influences  the  rms  normal  and  longitudinal  velocities.  Both  figures 
indicate  that  the  amplitude  increases  as  the  frequency  increases.  In 
these  figures,  the  frequency  is  increased  by  a  factor  of  16  from 
to*  1/8  to  to*  2.  The  points  of  minimum  amplitude  also  shift  upward  as 
the  frequency  increases.  It  may  be  useful  to  normalize  in  other  ways 
to  reduce  this  dependence1  of  the  amplitude  in  the  freestream  on 
frequency. 

A  possible  explanation  for  the  larger  amplitudes  as  frequency 
increases  is  based  on  the  analytical  solution  (A.  9)  of  Appendix  A  for 
the  viscous,  uniform  flow  case.  The  amplitude  is 


(A. 9) 
(A. 8) 


Hence,  the  coefficient  of  sinay  is  proportional  to  the  square  root  of 
frequency  when  p /&  «<1.  rFor  other  cases,  eqn.(A.8)  above  also 
suggests  that  the  coefficient  to/B  would  increase  as  cU  increases.  We 
offer  this  explanation  as  tentative,  however,  because  tne  above 
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Figure  3.20  The  mean  velocity  profile  and  first  two  derivatives  for  a 
favorable  pressure  gradient  boundary  layer  with  a 
Falkner-Skan  parameter  0.5 
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Figure  3.21  The  variation  of  the  normal  velocity  fluctuation  with 
distance  from  the  wall  for  a  decaying  wave,  exp(-£x) ,  for 
Z3  =<^=  0.5  and  R^=1000  for  a  favorable  pressure  gradient 
boundary  layer  with  a  Falkner-Skan  parameter  0.5 
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Figure  3.22  The  variation  of  the  longitudinal  velocity  fluctuation 
with  distance  from  the  wall  for  a  decaying  wave, 
exp(-^x),  for  £  =  0.5  and  R^=1000  for  a  favorable 

pressure  gradient  boundary  layer  with  a  Falkner-Skan 
parameter  0.5 
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Figure  3.23  The  variation  of  the  vorticity  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave,  exp(-sx) r  for 
^  —  co  —  0.5  and  R$=1Q00  for  a  favorable  pressure  gradient 
boundary  layer  with  a  Falkner-Skan  parameter  0.5 
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Figure  3.24  The  variation  of  the  pressure  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave,  expC-Sx)  ,  for 
0.5  and  Rj=1000  for  a  favorable  pressure  gradient 
boundary  layer  with  a  Falkner-Skan  parameter  0.5 
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Figure  3.25  The  variation  of  the  Reynolds  stress  with  distance  from 
the  wall  for  a  decaying  wave#  exp(-^ax)  ,  for  (3  maJ  =  0.5 
and  R$  =1000  for  a  favorable  pressure  gradient  boundary 
layer  with  a  Falkner-Skan  parameter  0.5 
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Figure  3.26  The  variation  of  the  energy  production  with  distance  from 
the  wall  for  a  decaying  wave,  exp(-,0x>  ,  for  p  =a>  =  0.5 
and  =1000  for  a  favorable  pressure  gradient  boundary 
layer  with  a  Falkner-Skan  parameter  0.5 
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Figure  3.28  The  variation  of  the  normal  velocity  fluctuation  with 
distance  from  the  wall  for  a  decaying  wave,  exp(-^x) ,  for 
/3  =a>=  0.5  and  R5=1000  for  an  unfavorable  pressure 
gradient  boundary  layer  with  a  Falkner-Skan  parameter 
-0.05 
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Figure  3.29  The  variation  of  the  longitudinal  velocity  fluctuation 
with  distance  from  the  wall  for  a  decaying  wave, 
exp(-/3x)  ,  for  p  =co  =  0.5  and  Rg=10Q0  for  an  unfavorable 
pressure  gradient  boundary  layer  with  a  Falkner-Skan 
parameter  -0.05 


AEDC-TH-83-4 


510  03  Z 


xo.o 


8*0 


6.0 

Distance  from 
wall,  y 

4.0 


2.0 


0.0 

-50,0  -25.0  00.0  25.0  50.0 

Complex  amplitude  of  the  vorticity  fluctuation,  Z(y) 


Figure  3.30  The  variation  of  the  vorticity  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave,  exp(-jgx)  ,  for 
p**co**  0.5  and  Rj=1000  for  an  unfavorable  pressure 
gradient  boundary  layer  with  a  Falkner-Skan  parameter 
-0.05 
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Figure  3,31  The  variation  of  the  pressure  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave#  expf-jSx)  ,  for 
p=eOsa  0.5  and  R^wlOOO  for  an  unfavorable  pressure 
gradient  boundary  layer  with  a  Falkner-Skan  parameter 
-0.05 
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Figure  3.33  The  variation  of  the  energy  production  with  distance  from 
the  wall  for  a  decaying  wave,  exp(-px)  ,  for  (3  **aj=  0.5 
and  Rg  =1000  for  an  unfavorable  pressure  gradient  boundary 
layer  with  a  Falkner-Skan  parameter  -0.05 
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Figure  3.34  Effects  of  frequency  on  the  rms  normal  velocity 
fluctuation,  with  the  other  parameters  held  constant. 
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Figure  3.35  Effects  of  frequency  on  the  rms  longitudinal  velocity 
fluctuationr  with  the  other  parameters  held  constant. 
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argument  seems  unable  to  provide  any  guidance  on  the  influence  of 
Reynolds  number,  as  discussed  later. 

The  solution  for  the  velocity  fluctuation  in  the  viscous  sublayer, 
derived  in  Appendix  B,  shows  that  the  ratio 

1A2 

sublayer  thickness/boundary  layer  thickness  =  1/  ) ' 

Hence,  as  the  frequency  is  increased,  the  shorter  times  available  for 
viscous  diffusion  to  take  place  during  one  oscillation  cycle  leads  to 
thinner  sublayers.  Calculations  of  the  rms  velocity  for  five 
frequencies  differing  by  a  factor  of  5  are  shown  in  Figure  3.36. 
These  calculations  show  the  thinning  of  the  sublayer  as  frequency 
increases. 


3.2.4  Effects  OL  Reynolds  Number 

Calculations  of  the  rms  velocities  with  Reynolds  numbers  varying 
by  a  factor  of  10  (from  R$  =  100  to  1000)  are  shown  in  Figures  3.37  and 
3.39.  Clearly,  the  velocities  are  rather  insensitive  to  the  Reynolds 
number,  except  in  the  sublayer. 

From  our  discussion  of  the  previous  section,  one  would  expect  that 
one  effect  of  increasing  the  Reynolds  number  would  be  to  thin  the 
sublayer ,  S*.  /  &  =  Evidence  that  this  is  true  is  shown  in 
Figure  3.39  where  the  rms  longitudinal  velocity  is  plotted  in  the 
sublayer  region  for  nine  Reynolds  numbers  differing  by  a  factor  of  20. 

In  Figure  3.39,  also  note  the  apparent  coalescence  of  the  curves 
above  the  sublayer,  which  further  supports  our  conclusion  that  the 
velocities  are  rather  insensitive  to  R^  except  in  the  sublayer. 


3.2.5  Effects  £f  Reynolds  Number  and  Frequency  When  a> Rc  la  Held 
Constant. 

Since  the  ratio  of  sublayer  thickness  to  boundary  layer  thickness 
is  ( [oRjT1/ ,  then  if  the  frequency  and  Reynolds  number  are  changed  but 
their  product  is  held  constant,  then  one  would  expect  the  sublayer  to 
be  relatively  unchanged.  To  test  this  hypothesis,  we  calculated  the 
following  cases 


=  2.0 
=»  1.0 
as  =  0.5 
=  0.25 
OJ  =  0.125 


Ri=  250 
Rj=  500 
Rfi=  1000 
R$=  2000 
R4=  4000 


u>  500 
«^RS=  500 
<oRl=  500 
«'Rs=  500 
^  R^=  500 


The  rms  longitudinal  velocity  in  the  sublayer  is  plotted  in  Figure  3.40 
for  these  five  cases.  The  sublayer  thickness  appears  to  be  unchanged 
although  the  frequency  and  Reynolds  number  separately  have  been  changed 
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Figure  3.36  Effect  of  frequency  on  the  rms  longitudinal  velocity 
fluctuation  in  the  sublayer,  with  the  other  parameters 
held  constant. 
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Figure  3.37  Effects  of  Reynolds  number  on  the  rms  normal  velocity 
fluctuation,  with  the  other  parameters  held  constant. 
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Figure  3.38  Effects  of  Reynolds  number  on  tne  rms  longitudinal 
velocity  fluctuation,  with  tne  other  parameters  held 
constant. 


AEDC-TH-83-4 


546  05  URMS 


Figure  3.39  Effect  of  Reynolds  number  on  the  rms  longitudinal 
velocity  fluctuation  in  the  suolayer.  Note  tne 
coalesence  of  the  curves  into  the  (nearly)  Reynolds 
number  independent  curve  outside  of  the  sublayer. 
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Figure  3.40  Effects  of  changing  both  the  frequency  and  the  Reynolds 
number,  but  holding  their  product  constant, 500. 
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by  a  factor  of  16. 


3.2.6  Effects  of  v-wavenumher  or  Decay  Bata 

To  obtain  evidence  of  the  influence  of  the  y-wavenumber  (or  decay 
rate  in  the  streamwise  direction) ,  we  calculated  the  rms  velocity  for 
(3=0.10,  0.25,  0.50,  and  0.75  with  the  other  parameters  fixed.  These 
results  are  plotted  in  Figure  3.41.  The  higher  frequency  oscillations 
associated  with  increasing  the  y-wavenumber  are  apparent.  Also,  the 
maximum  amplitude  increases  as  jS  decreases.  Based  upon  the  guidelines 
of  Appendix  A,  the  coefficient  of,  sin fly  is 

'W/z 

Hence,  as  p  decreases,  this  coefficient  increases,  in  agreement  with 
our  calculations  for  a  Blasius  layer. 


3.3  NUMERICAL  RESULTS  FOR  GROWING  STANDING  WAVES 


3.3.1  Inviscid  Relationship  Between  the  Solutions  for  BfiCayinq  and 
Growing  Waves 

The  Orr-Sommerfeld  equation  (eqn.  2.6b)  cam  be  separated  into  the 
following  two  real  equations: 

(1)  Equation  originating  from  the  real  part: 

{  3(&f)  -  -  fa  (r y  *  =  «.13, 

(2)  Equation  originating  from  the  imaginary  part: 

where  r=+l  or  -1  for  growing  or  decaying  waves,  respectively. 

If  ^(y)  ■  #(y)+i<5f-(y>  is  a  solution  for  a  decaying  wave  (r=-lJ, 
then  we  ask:  Is  the  complex  conjugate  solution  <?  =#-i a  solution  for 
the  growing  wave  case  ( r*>+l) ?  If  one  substitutes  r=-l  with  <P  »(£+i$, 
and  if  one  substitutes  r»+l  with  ,  one  finds  the  two  systems 

are  identical  only  i£  the  flows  are  inviscid.  The  viscous  term  does 
not  transform  in  that  manner. 

As  we  shall  see  from  numerical  calculations  for  growing  and 
decaying  waves,  this  inviscid  rule,  well  known  from  inviscid  stability 
theory,  is  a  practical  guideline  for  the  high  Reynolds  number  viscous 
case  as  well. 
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3.41  Effect  of  the  y-wavenumber  or  decay  rate  on  the  rms 
normal  velocity  fluctuation,  with  tne  other  parameters 
held  constant. 
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3.3.2  Bias ins  case  with  growing  waves 

Figures  3.42-3.47  present  numerical  results  for  a  wave  growing 
exponentially  in  the  streamwise  direction  as  exp(+^x) .  These  figures 
can  be  compared  with  the  results  for  the  decaying  wave  as  indicated 
below. 

exp(-px)  wave  exp(+px)  wave 


<P  Figure  3.2  3.42 

f  3.3  3.43 

Z  3.6  3.44 

77  3.11  3.45 

-uv  _  3.17  3.46 

-uvdU/dy  3.18  3.47 


While  <P  for  the  exp(+/Sx)  wave  is  not  exactly  the  complex 
conjugate  of  <P  for  the  exp(-^x)  wave,  the  inviscid  theorem  appears  to 
be  a  reasonable  guideline  for  the  large  Reynolds  number  cases.  The 
longitudinal  amplitude,  f=fr+if(-,  for  the  decaying  wave  is 
approximately  -fr+if/  for  the  growing  wave  outside  of  the  sublayer. 
The  overshoot  of  fr  at  the  sublayer  edge  for  a  decaying  wave  appears  to 
be  similar  to  the  overshoot  of  t;  for  the  growing  wave.  A  significant 
discrepancy  from  the  inviscid  guideline  arises  for  the  pressure. 

The  Reynolds  stress  in  Figure  3.46  has  a  reversal  of  sign  compared 
with  that  stress  in  Figure  3.17.  This  sign  change  leads  to  a  sign 
change  in  the  energy  production.  Figure  3.47  as  compared  with  3.18.  A 
positive  energy  production  indicates  that  energy  is  being  transferred 
from  the  mean  flow  to  the  disturbance  flow. 


3.4  NUMERICAL  RESULTS  WHEN  THE  NO-SLIP  CONDITION  AT  THE  WALL  IS 
REPLACED  BY  A  NO-SHEAR  CONDITION. 

To  give  some  indication  of  the  role  of  the  viscous  sublayer,  in 
Figures  3.48-3.51  we  replaced  the  no-slip  condition,  u  =  0  or  <£=  0, 
with  a  no-shear  condition  at  the  wall  9 


=  O  O'”  ~  °  (3.15) 

These  figures  with  a  no-shear  condition  may  be  compared  with  the 
corresponding  ones  with  a  no-slip  condition  as  tabulated  below: 


no-slip 

no-shear 

Figure  3.2 

3.48 

f 

3.3 

3.49 

z 

3.6 

3.50 

77 

3.11 

3.51 
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Complex  amplitude  of  the  normal  velocity  fluctuation,  ?(y) 


Figure  3.42  The  variation  of  the  normal  velocity  fluctuation  witn 
distance  from  the  wall  for  a  growing  wave,  exp(+^x) , 
for  1 3-<o=  0.5  and  Rg=1000  with  a  Blasius  boundary  layer. 
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Figure  3.43  The  variation  of  the  longitudinal  velocity  fluctuation 
with  distance  from  the  wall  for  a  growing  wave,  exp(+px) , 
with  0.5  and  R$=1000  with  a  Blasius  boundary 

layer. 
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Figure  3.45  The  variation  of  the  pressure  fluctuation  with  distance 
from  the  wall  for  a  growing  wave,  exp(+px) ,  for 
0.5  and  Rj=1000  with  a  Blasius  boundary  layer. 
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Figure  3.46  The  variation  of  the  Reynolds  stress  with  distance  from 
the  wall  for  a  growing  wave,  exp(+px)  ,  for  0.5  and 

R^=*1000  with  a  Blasius  boundary  layer.  * 
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Figure  3.47  The  variation  of  the  energy  production  with  distance  from 
the  wall  for  a  growing  wave,  exp(+|Sx)  ,  for  s  = a. 0.5  and 
R£=1000  with  a  Blasius  boundary  layer. 
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Figure  3.48  The  variation  of  the  normal  velocity  fluctuation  witn 
distance  from  the  wall  for  a  decaying  wave,  exp(-^x) ,  for 
|3  =cu=  0.5  and  R$=100Q  for  a  Blasius  layer,  and  where  the 
no-slip  condition  on  the  fluctuating  flow  has  been 
replaced  by  a  no-shear  condition. 
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The  impermeability  condition  at  the  wall  was  imposed  for  both  sets  of 
figures. 

As  expected  and  shown  clearly  in  Figure  3.49  for  the  longitudinal 
velocity,  no  sublayer  forms  when  the  no-slip  condition  is  replaced  by  a 
no-shear  condition.  For  the  variables  plotted,  the  values  above  the 
sublayer  region  in  Figures  3.2,  3.3,  3.6,  and  3.11  are  nearly  the  same 
as  the  values  in  Figures  3.48-3.51. 

We  conclude  that  the  sublayer  is  relatively  passive,  i.e.,  it  is 
driven  by  the  external  flow  above  the  sublayer,  and  that  external  flow 
is  not  influenced  very  greatly  by  the  sublayer.  Except  in  the  sublayer 
region  itself  of  thickness  0(1  /  ,  the  results  indicate  that  if  the 
Reynolds  number  and  frequency  are  large,  »i,  then  the 
replacement  of  the  no-slip  condition  by  a  no-shear  condition  has  a 
relatively  small  influence  on  the  solutions.  This  result  differs 
sharply  from  the  active  role  which  the  sublayer  takes  for  the 
Tollmien-Schlichting  wave. 


82 


10.0 

6.0 

6.0 

Distance  from 
wall,  y 

4.0 

2.0 

0.0 

-2.5  -2.0  -1.5  -1.0  -0.5  0.0  0.5  1.0  1.5  2.0  2.5 

Complex  amplitude  of  the  longitudinal  velocity  fluctuation,  f(y) 


519  02  F 


Figure  3.49  The  variation  of  the  longitudinal  velocity  fluctuation 
with  distance  from  the  wall  for  a  decaying  wave, 
exp{-0x)  ,  for  (3  ~a/s  0.5  and  R^=1000  for  a  Blasius  layer, 
and  where  the  no-slip  condition  on  tne  fluctuating  flow 
has  been  replaced  by  a  no-shear  condition. 
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Figure  3.50  The  variation  of  the  vorticity  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave,  exp(-px) ,  for 
puetJm  0.5  and  R^=1000  for  a  Blasius  layer,  and  where  the 
no-slip  condition  on  the  fluctuating  flow  has  been 
replaced  by  a  no-shear  condition. 
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Figure  3.51  The  variation  of  the  pressure  fluctuation  with  distance 
from  the  wall  for  a  decaying  wave,  exp(-px),  for 
p  =&=  0.5  and  R^=<1000  for  a  Blasius  layer,  and  where  the 
no-slip  condition  on  the  fluctuating  flow  has  been 
replaced  by  a  no-shear  condition. 
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4.0  SUMMARY,  DISCUSSION  AND  CONCLUSIONS 


4.1  INTRODUCTORY  COMMENTS 

This  report  documents  the  analytical  and  numerical  evidence  of  the 
existence  of  a  pair  of  standing  waves  in  a  boundary  layer.  These  waves 
oscillate  sinusoidally  in  time,  but  because  of  the  purely  imaginary 
x-wavenumber  and  phase  speed,  the  waves  do  not  travel.  The  amplitudes 
vary  exponentially  in  the  streamwise  direction  without  any  spatial 
oscillation.  Both  growing  and  decaying  solutions  are  possible.  Far 
above  the  boundary  layer,  these  irrotational  fluctuations  oscillate 
neutrally  in  the  y-direction. 

Some  evidence  of  the  existence  of  these  waves  is  based  on  several 
numerical  solutions  of  the  Orr-Sommerfeld  equation,  analytical 
solutions  of  simplified  forms  of  the  Orr-Sommerfeld  equation,  and  an 
asymptotic  solution  of  the  Orr-Sommerfeld  equation.  The .asymptotic 
solution  by  Shunichi  Tsuge  as  applied  to  the  stability  problem  is 
summarized  in  Ref. 26.  The  use  of  that  asymptotic  solution  in  the 
spatial  initial-value  problem  will  be  summarized  in  another  report. 

The  Orr-Sommerfeld  equation  describes  the  evolution  of 
small-amplitude  fluctuations  in  a  boundary  layer  growing  sufficiently 
slowly  such  that  it  is  meaningful  to  represent  it  as  a  layer  of  uniform 
thickness.  The  characteristics  of  the  boundary  layer  enter  into  the 
equation  as  two  coefficients.  The  Falkner-Skan  family  of  boundary 
layers  was  used  so  that  evidence  could  be  obtained  for  fluctuations  in 
boundary  layers  with  favorable,  zero,  and  unfavorable  mean  pressure 
gradients.  Also,  since  roughness  can  influence  the  mean  profile,  these 
profiles  provide  some  indication  of  how  roughness  can  influence  the 
standing  waves.  Other  profiles  will  be  used  later. 

To  illustrate  how  standing  waves  originate,  the  interaction 
between  freestream  vortical  fluctuations  and  leading  and  trailing  edges 
of  a  flat  plate  are  summarized  in  Chapter  1.  Those  studies,  while 
including  the  features  of  leading  and  trailing  edges,  do  not  include 
all  of  the  features  of  the  Orr-Sommerfeld  equation.  However,  the 
problems  share  three  essential  features  of  the  fluctuations  in  tne 
freestream:  (1)  they  oscillate  in  time,  (2)  decay  or  grow  in  the 
streamwise  direction,  and  (3>  are  irrotational.  This  report  documents 
the  evidence  that  the  Orr-Sommerfeld  admits  solutions  with  those 
characteristics,  and  that  other  features  arise  because  of  the  shearing 
mean  flow  and  viscosity.  The  problem  of  predicting  the  origin  of  the 
various  waves  is  a  topic  of  continuing  research. 


4.2  SUMMARY  OF  THE  FORMULATION  AND  NUMERICAL  TECHNIQUE 

For  small-amplitude  disturbances  in  shear  layers  which  develop 
sufficiently  slowly  such  that  the  parallel-flow  representation  is 
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adequate,  the  objective  is  to  obtain  solutions  of  the  equation 

{h*  DMTx  *V*}  7V‘  =  O  u.i) 

which  are  of  the  Fourier-Laplace  form 

v(x,yj)=  <P(y)e  ?  (4.2, 

where  ^  and  uf  are  pure  real.  These  standing-wave  solutions  are 
exponentially  growing  or  decaying  in  the  streamwise  direction,  and 
oscillate  sinusoidally  in  time.  The  motivations  for  seeking  solutions 
of  this  form  are  that  such  waves  appear  in  several  analytical  solutions 
where  the  mean  flow  is  uniform.  Furthermore,  in  an  asymptotic  analysis 
of  an  initial-value  problem  with  a  boundary  layer,  two  mathematical 
poles  indicate  that  waves  of  this  form  would  exist  in  boundary  layers. 
For  solutions  of  this  form,  the  x-wavenumber  is 

o (  =  -t'r/S  (4.3) 

m  =  +1  for  a  growing  wave 
r  =  -1  for  a  decaying  wave 


The  complex  phase  speed,  which  is  pure  imaginary,  is 


(4.4) 


Letting  oi  =-ir ,  the  Orr-Sommerfeld  equation  is 

{  (fm- 

The  impermeability  and  no-slip  boundary  conditions  at  the  wall, 
and  the  boundedness  condition  of  the  solution  far-away,  are 

-  O  ctf  y~0  *t  <jP  es  as y  —■  <**  (4.5b,c,d) 


When  0=1  and  (L^=  0,  the  Orr-Sommerfeld  equation  has  the 
solutions  " 


<?(</)  = 

A  e~ 

.  /7py 

■+ 

&e  +  C  cos^Sy 

*  Ossify 

(4 

.€) 

or  equivalently 

<?(</)  ~ 

Ae~ 

(4. 

.7) 

where  the  exponent  is 

m  - 

(rP*rP 

(4. 

.8) 

For  the  boundary 

layer  case,  B  =  0. 

For  most  cases,  we  normalize 

by 

setting  C=1  in  eqn.(4.6).  In  retrospect,  the  analytical  results  of 
Appendix  A  and  the  calculations  suggest  that  setting  the  coefficient 
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D=1  would  be  a  better  choice,  and  might  lead  to  amplitudes  of  the 
freestream  disturbances  to  be  less  sensitive  to  the  parameters. 
However#  for  any  normal  modes  analysis,  the  amplitude  ^(y)  can  be 
multiplied  by  a  complex  constant  and  rescaled.  The  essential 
requirement  in  carrying  out  calculations  for  the  initial-value  problem 
is  consistency  between  the  normalizations  of  the  normal  modes 
calculations  and  the  initial-value  theory. 

Pa r-away  from  the  boundary  layer  edge  where  the  term  Aexp(-my)  is 
negligible,  the  flowfield  varies  sinusoidally 

,  ,  /  i  ±  fix- reef 

=  (  (osfy  PySnfyJe  \  (4>9) 

This  freestream  disturbance  ia  lr  rotational.  The  constant  D  is 
obtained  numerically. 

In  general,  the  constant  D  and  the  solution  <^(y)  of  trie 
Orr-Soramerfeld  equation  depend  on  the  frequency,  y-wavenuraber,  Reynolds 
number,  mean  velocity  profile,  and  the  choice  of  whether  the  wave  is 
growing  or  decaying.  The  effects  of  the  mean  pressure  gradient  and  one 
of  the  effects  of  surface  roughness  can  be  introduced  tnrough  the  mean 
velocity  profile,  UCy>. 

Numerical  solutions  of  the  Orr-Sommerfeld  equation  are  obtained  by 
an  expansion  in  a  series  of  (typically  48)  Chebyshev  polynomials.  The 
coefficients  of  those  polynomials  are  obtained  as  solutions  of  a  set  of 
linear,  algebraic  equations.  The  matrix  of  coefficients  is  reduced  by 
Gauss-Jordan  elimination. 

Calculations  are  carried  out  on  a  Digital  Equipment  Corporation 
PDP-11/23  computer  using  64K  bytes  of  memory  and  another  64K  bytes  of 
extended  memory  for  virtual  arrays.  The  disk-over layed  FORTRAN  .IV 
program  executes  in  about  145  seconds  when  using  the  KEF11— AA  floating 
point  chip,  beginning  with  data  for  the  mean  velocity  profile,  and 
ending  with  the  coefficients  of  the  polynomials  and  the  amplitude  ^(y) • 
The  program  executes  in  about  80  seconds  when  using  the  FPF-11  floating 
point  processor. 


4.3  SUMMARY  OF  RESULTS,  DISCUSSION  AND  CONCLUSIONS 

The  calculations  for  the  standing  waves  indicate  that  the 
longitudinal  velocity  fluctuation  can  bn  about  as  large  inside  the 
boundary  layer  as  the  fluctuation  in  the  freestream.  This  conclusion 
differs  significantly  from  the  solutions  with  decaying  vortices  in  the 
freestream,  as  can  be  seen  for  the  rms  longitudinal  velocity  in 
Figure  4.1  from  Ref.  15.  This  plot  has  the  same  frequency, 
y-wavenumber,  Reynolds  number  and  mean  velocity  as  for  Figure  3.14,  but 
the  velocities  deep  inside  the  boundary  layer  are  small  in  Figure  4.1. 
The  corresponding  plot  for  rms  vorticity  from  Ref. 15  is  shown  in  Figure 
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Figure  4.1  The  rms  longitudinal  velocity  for  the  case  of  spatially- 
decaying  vortical  fluctuations  for  R$  =  1000,  p  =  co~  0.5, 
and  a  Blasius  boundary  layer.  Compare  this  figure  for 
a  travelling  wave  with  Figure  3.14  for  a  standing  wave. 
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4.2.  This  plot  shows  that  a  layer  of  vorticity  exists  near  the 
boundary  layer  edge.  The  phase  of  this  vorticity  is  such  that  the  flow 
induced  by  this  layer  of  vorticity  tends  to  cancel  the  flow  induced  by 
the  freestream  vortical  fluctuations.  This  cancellation  process  does 
not  appear  to  take  place  with  the  standing  waves.  Inviscid,  analytical 
solutions  of  the  Rayleigh  equation  with  broken-line  mean  velocity 
profiles  would  help  clarify  this  process.  However ,  tne  author  has 
learned  from  several  studies  that  tne  2-D  response  in  boundary  layers 
is  sensitive  to  the  phase  speed  of  freestream  disturbances. 

The  rms  longitudinal  velocity  for  a  Tollmien-Schlichting  wave  is 
plotted  in  Figure  4.3  for  comparison  with  the  other  waves. 

Although  the  "freestream  disturbance"  for  tne  standing  wave  is 
ir rotational,  the  fluctuating  flow  inside  the  boundary  layer  is 
rotational.  That  vorticity  originates  primarily  from  two  sources:  the 
production  of  fluctuating  vorticity  inside  the  boundary  layer  by  the 
term  -vUyy  ,  and  by  diffusion  in  the  y-direction,  £  ,  to  and  from  the 
wall  in  tne  unsteady  viscous  sublayer. 

The  main  features  of  a  thin  viscous  sublayer  can  be  described  by 
the  unsteady,  forced  diffusion  equation  (B.l) .  The  forcing  arises  from 
the  pressure  gradient  impressed  on  the  sublayer  by  the  flow  above  the 
sublayer. 

If  the  no-slip  condition,  u  »  0  or  <Pa  =  0,  is  replaced  by  a 
no-shear  condition  at  the  wall,  uy  =  0  or  =  0,  calculations  show 
that  the  disturbance  is  not  significantly  affected  (except  for  the 
sublayer  itself) .  The  sublayer  does  not  form  when  tne  no-shear 
condition  is  imposed.  The  result  that  the  solution  outside  the 
sublayer  region  is  relatively  unchanged  can  be  interpreted  as  the 
sublayer  being  relatively  passive  under  the  high  Reynolds  number ,  high 
frequency  condition. 

A  remarkable  balance  takes  place  for  these  standing  waves.  The 
fluctuating  vorticity  decays  (or  grows)  exponentially  in  the  streamwise 
direction  and  oscillates  in  time  in  exactly  the  same  fashion  as  tne 
irrotational  fluctuations  in  the  freestream  disturbance.  The  vorticity 
induces  a  flowfield  which  generates  additional  vorticity  through  the 
production  term,  and  that  flowfield  also  affects  the  diffusion 
processes  in  the  sublayer.  Also,  the  impermeability  condition 
essentially  sets  up  an  image  system.  The  "history"  of  the 
flowfield,  including  the  oscillations  and  phasing,  is  important  since 
the  vorticity  convecting  downstream  in  the  shear  flow  must  arrive  at 
each  station  and  be  altered  by  the  production  and  diffusion  in  a  such  a 
manner  that  the  vorticity  field  is  also  a  standing  wave.  With  these 
varied  processes  taking  place  simultaneously,  it  is  amazing  that  a 
particular  arrangement  is  possible  such  that  the  flowfield  varies 
exponentially  in  space  and  oscillates  in  time. 

Not  all  irrotational  fluctuations,  which  decay  or  grow 
exponentially  in  the  freestream  and  oscillate  in  time,  are  solutions  of 
the  Orr-Sommerfeld  equation,  although  they  may  be  solutions  of  the 
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Figure  4.2  The  rms  fluctuating  vorticity  for  the  case  ot  spatially- 
decaying  vortical  fluctuations  for  R4  =  1000,  /3-^=  0.5, 
and  a  Blasius  boundary  layer. Note  the  bulge  in  vorticity 
near  the  boundary  layer  edge.  The  bulge  further  out  is 
representative  of  the  vortical  fluctuations  in  the 
f reestream. 
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Figure  4.3  The  rms  longitudinal  velocity  for  a  Tollmien-Schlichting 
wave  for  R$  =  821.835,  <*=  (0.2888,  0.0),  and  com  0.1074 
with  a  Blasius  boundary  layer. 
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partial-differential  equation  (2.1).  in  the  boundary  layer,  the  flow 
may  be  much  more  complicated  then  the  assumed  solution  form  (2. 3, 2. 4). 

An  important  idea  from  the  initial -value  theories,  however,  is 
that  the  more  complicated  solutions  can  be  constructed  by  superposition 
of  the  basic  Fourier-Laplace  solutions.  In  some  cases,  tne 
superposition  is  a  summation  of  waves  of  finite  amplitude.  In  other 
cases,  the  superposition  is  an  integral  over  waves  of  infinitesimal 
strengths . 


The  other  Fourier-Laplace  solutions,  including  the  stability  waves 
and  the  two  spatial-solutions  with  vortical  fluctuations  in  tne 
freestream,  also  have  remarkable  balances  which  lead  to  solutions  with 
simple  functional  dependences  on  x  and  t. 


From  the  theory  of  inviscid  stability  waves,  it  is  known  that  if  <p 
is  the  solution  for  the  wavenumber  «*',  then  ^'is  tne  solution  for  the 
complex  conjugate  wavenumber , W *.  Calculations  based  on  the  viscous 
Orr-Sommerfeld  equation,  however,  indicate  that  tne  inviscid 
relationship  is  a  reasonable  guideline  for  relating  the  solutions  for 
the  viscous  decaying  and  growing  standing  waves.  The  success  is 
probably  related  to  the  passive  character  of  the  sublayer  when  the 
sublayer  is  thin,  and  the  general  insensitivity  of  the  standing  wave 
solutions  to  the  Reynolds  number  (except  in  the  sublayer) .  The  details 
of  the  flow  in  the  sublayer  do  not  transform  according  to  inviscid 
guidelines.  The  wavenumber  for  the  growing-wave  case,  «*<  -  -ia  is  the 
complex  conjugate  of  that  wavenumber  for  the  decaying  case,**'  ■  +ip  . 

The  functional  dependence  of  the  standing  waves  indicates  that  the 
waves  do  not  travel  in  the  streamwise  direction*  However,  a  movie 
shows  that  wave  motion  takes  place  in  the  y-direction.  In  the  bounds ry 
layer,  that  motion  is  not  a  simple  sinusoidal  wave,  and  it  does  not 
travel  in  the  y-direction  with  a  constant  speed.  Far  above  tne 
boundary  layer  where  the  only  remaining  influence  of  the  boundary  layer 
is  through  the  value  of  the  constant  D,  eqn.(4.9)  can  be  rearranged  as 


This  expression  indicates  that  the  irrotational  freestream  disturbance 
is  a  superposition  of  outgoing  and  incoming  waves.  The  amplitudes  of 
these  waves  are  different  since,  in  general,  E  ^  1. 

All  of  the  available  information  is  consistent  with  the  conclusion 
that  the  two  unsteady,  exponentially-varying  standing  waves  are  basic, 
independent,  Fourier-Laplace  solutions  of  the  viscous,  parallel-flow 
Orr-Sommerfeld  equation  describing  small-amplitude  fluctuations  in  a 
boundary  layer.  This  evidence  is  based  on  several  analytical 
solutions,  a  formulation  of  the  spatial  initial-value  problem  using  an 
asymptotic  solution  of  the  Orr-Sommerfeld  equation,  and  the  numerical 
calculations  with  adverse,  zero,  and  favorable  pressure  gradient 
boundary  layers.  For  mathematical  completeness,  these  two  solutions 
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join  the  list  consisting  of  the  discrete  set  of  stability  waves  and  the 
pair  of  solutions  accounting  for  the  propagation  of  vortical 
fluctuations  in  the  freestream.  While  this  report  has  focused  on  the 
standing  waves  in  boundary  layers ,  the  author  believes  that  equivalent 
waves  arise  in  free-shear  layers ,  jets,  and  wakes. 


4.4  DESCRIPTION  OF  PROBLEMS  WITH  IRROTATIONAL  DISTURBANCES  IN  THE 
FREESTREAM 

An  objective  of  this  report  is  to  document  the  evidence  that  a 
pair  of  standing  wave  solutions  exists  which  have  the  characteristics 
of  i rrotattonal  fluctuations  "far-away"  fiom  the  boundary  layer .  In 
all  results,  the  vorticity  and  velocities  satisfy  parallel-flow 
equations.  Those  solutions  summarized  in  Chapter  1  and  in  Appendix  A 
all  have  a  uniform  mean  flow.  Consequently,  there  can  be  no  Rayleigh, 
Tollmien-Schlichting,  or  other  instability  waves  in  those  problems. 
Yet,  both  rotational  and  irrotational  disturbances  are  present.  The 
irrotational  solutions  fall  into  two  classes: 

(1)  Travelling  waves  which  arise  through  the  impermeability 
condition  of  the  plate.  In  more  sophisticated,  viscous 
solutions  with  shearing  flows,  these  travelling  waves  are  a 
part  of  the  response  of  a  boundary  layer  to  vortical 
fluctuations  in  the  freestream.  This  irrotational  portion 
does  not  survive  far-away  from  the  boundary  layer. 

(2)  Standing  waves,  which  can  also  be  set  up  by  the 
impermeability  condition,  but  which  do  not  travel.  Unlike 
the  vortical  fluctuations  which  have  no  pressure  fluctuations 
in  the  freestream,  these  standing  waves  have  (sometimes 
strong)  pressure  fluctuations  which  can  be  estimated  from 
linear  equations.  This  flowfield  is  not  describable  by 
stability  waves,  and  is  not  describable  by  the  solutions  (1) 
above. 

An  initial-value  problem  formally  illustrates  that  the  two 
standing  waves  are  entities  distinct  from  the  poles  of  the  stability 
waves  and  from  the  branch  lines  of  the  vortical  fluctuations.  The 
standing  waves  in  this  report  are  described  by  tne  viscous 
Or r- Sommer f eld  equation  with  a  boundary  layer,  rather  than  the  uniform 
flows  as  used  in  Chapter  1  and  Appendix  A.  The  formal  linking  of 

(a)  the  standing  waves  set  up  by  the  interactions  of  vortical 
fluctuations  with  the  leading  edge,  as  analyzed  in  Refs. 6  and  7,  and 

(b)  the  standing  waves  described  by  the  Orr-Somraerfeld  equation  is  a 
topic  of  continuing  research  which  is  not  likely  to  be  resolved 
rigorously  in  the  near  future.  The  information  presently  available  is 
based  on  a  sequence  of  problems:  (a)  solutions  of  the  Orr-Sommerfeld 
equation,  (b)  uniform-flow,  viscous  solutions,  (c)  uniform-flow 
inviscid  solutions,  and  (d)  uniform-flow  inviscid  solutions  with  a 
leading  edge.  In  each  of  the  problems  (a,b,c),  the  standing  waves 
reappear  in  ever-simplier  forms.  The  two  problems  (c,d)  have  identical 
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equations  and  boundary  conditions  downstream  of  the  leading  edge,  and 
both  problems  have  standing  waves.  In  tne  future,  the  solution  from 
problem  (d)  will  be  Fourier-analyzed  to  show  how  the  solutions  from 
problem  (c)  formally  enter  into  problem  (d) .  The  author  has  little 
doubt  that  such  a  link  is  possible  and  will  serve  as  a  guideline  to 
analyzing  the  more  complicated  viscous  cases. 

The  presence  of  irrotational  fluctuations  around  airfoils  and  in 
turbomachinery  is  no  surprise  to  unsteady  aerodynamicists.  Potential 
fluctuations  appear  in  the  unsteady  aerodynamics  theories  of  airfoils 
encountering  gusts,  airfoils  oscillating  in  angle  of  attack,  airfoils 
oscillating  in  the  direction  perpendicular  to  the  freestream,  and 
unsteady  flap  motions.  Experimental  data  confirms  the  presence  of  such 
fluctuations.  Potential  fluctuations  arise  when  rotor  blades  pass  by 
stator  blades  in  compressors  and  turbines.  While  some  of  these 
irrotational  fluctuations  are  travel ling -waves,  tne  unsteady 
circulation  about  an  airfoil  is  a  standing  wave  for  that  airfoil.  The 
flows  induced  by  the  vortex  sheet  shed  from  the  trailing  edge  are 
irrotational  (except  at  the  sheet  itself) .  These  flows  also  have 
standing  and  travelling  wave  components,  where  the  standing  wave  arises 
because  of  the  interaction  between  the  vortex  sheet  and  the  airfoil 
itself.  Far-downstream  of  the  airfoil,  the  principal  contribution  is 
the  travelling  wave  induced  by  the  vortex  sheet  itself. 

The  model  of  an  airfoil  oscillating  in  a  direction  perpendicular 
to  the  flow  suggests  that  some  of  the  linear  features  model 
vibration  at  low  frequencies  would  appear  as  standing  waves.  The 
oscillating  dipole  character  of  an  oscillating  cylinder  beside  a 
boundary  layer  is  a  unsteady  standing  wave.  Oscillating  airfoils, 
oscillating  source-sinks,  and  airfoils  with  oscillating  flaps  which  are 
positioned  in  the  vicinity  of  a  plate  with  a  boundary  layer  will  also 
produce  standing  waves.  Travelling  waves  of  various  forms  will  also 
arise  in  these  cases. 

The  standing  waves  appear  to  be  one  of  the  downstream  influences 
of  the  leading  edge,  and  one  of  the  upstream  influences  of  a  trailing 
edge  of  a  body  interacting  with  freestream  vorticity  fluctuations.  The 
diffusive  wave  which  travels  upstream,  as  described  in  Refs.  17  and  27, 
is  another  mechanism  for  upstream  influence  in  this  elliptic  problem. 
Rockwell  (Ref.  25)  comprehensively  reviews  the  complex  and  interesting 
unsteady  flows  arising  when  vortical  fluctuations  interact  with  leading 
edges. 


Refs.  28-35  are  examples  of  numerical  solutions  which  have  been 
obtained  for  unsteady  boundary  layer  problems  with  finite  difference, 
vortex  filament,  and  spectral  techniques.  These  investigators  have 
been  faced  with  the  difficulty  of  properly  posing  the  boundary 
conditions  on  the  downstream,  outer,  and  upstream  boundaries  of  their 
calculational  regions.  This  problem  with  boundary  conditions  will 
become  more  severe  as  the  disturbances  adopt  the  characteristics  ot 
freestream  turbulence  and  its  interaction  with  tne  leading  edge,  as  the 
amplitudes  increase,  and  as  3-D  and  other  significant  features  are 
included. 
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The  standing  waves  could  be  an  important  factor  in  obtaining 
numerical  solutions  in  a  rectangular  domain,  since  those  waves  are  a 
mechanism  by  which  the  upstream  and  downstream  boundary  conditions  can 
influence  the  flow  elliptically.  Improperly  formulated  boundary 
conditions  could  excite  extraneous  fluctuations  in  the  calculational 
domain.  The  magnitude  of  these  unwanted  disturbances  would  depend  on 
the  amplitudes  at  the  boundaries,  as  well  as  the  y-wavenumber  which 
dictates  the  growth/decay  rate  in  the  domain,  and  the  length  of  the 
domain  itself.  The  reader  is  reminded  that  y-wavenumber  «  0  for 
non-decaying  waves  is  a  possible  limiting  case  of  these  standing  waves. 
Standing  waves  with  small  y-wavenumbers  would  influence  the  entire 
calculational  domain.  The  "purification"  of  a  solution,  by 
identification  and  elimination  of  unwanted  waves  in  the  calculation, 
could  be  a  useful  numerical  procedure. 

In  other  cases  as  suggested  by  the  paragraphs  above,  the  proper 
boundary  conditions  should  incorporate  standing  waves  which  naturally 
writer  into  and  leave  the  calculational  region. 

Possible  topics  for  future  research  include  the  investigations  of 
Cl>  3-D  standing  waves,  (2)  the  compressible  counterparts  to  tne 
present  incompressible  waves,  (3)  standing  waves  in  parallel-flow  jets, 
wakes,  and  free-shear  layers,  (4)  standing  waves  in  nonparallel 
shearing  layers,  and  (5)  standing  waves  in  boundary  layers  as  described 
by  the  unsteady  Prandtl  boundary  layer  equation.  More  detailed  studies 
are  needed  to  (6)  determine  the  influence  of  surface  roughness  on  these 
waves.  Further  development  of  (7)  the  initial-value  problem  would 
permit  an  initial  disturbance  to  be  decomposed  into  the  various 
standing  and  travelling  waves.  (8)  The  processes  which  create 
ir rotational  fluctuations  in  the  freestream  deserve  detailed  study. 
These  processes  include  the  interaction  of  vortical  freestream 
disturbances  with  leading  and  trailing  edges,  the  unsteady  motions  of 
airfoils  and  vibrations  of  models,  and  the  combined 
ir rotational/ rotational  disturbances  set  up  in  turbomachinery. 

In  closing,  many  additional  examples  can  be  cited  where  potential 
fluctuations  arise  in  practical  engineering  problems.  This  research  is 
intended  to  show  that  they  also  appear  as  solutions  of  tne 
Orr-Soramerfeld  equation  if  the  proper  parameters  and  exterior  boundary 
conditions  are  selected. 


96 


AEDC-TR-83-4 


REFERENCES 

1.  Rogler,  Harold  L.  and  Reshotko,  Eli.,  "Disturbances  in  a  boundary 

layer  introduced  by  vortex  array  representations  of  free  stream 
turbulence",  Workshop  on  Low  Speed  Boundary  Layer  Transition, 
Rand  Corporation,  Santa  Monica,  California  (16-17  July  1974); 
Abstract  appears  in  R-1752-ARPA/ONR  by  William  S.  King  (June 
1975) . 

2.  -  "Spatially  decaying  array  of  vortices".  The  Physics  of  Fluids, 

vol.  19,  no.  12  (December  1976). 

3.  - "Disturbances  in  a  boundary  layer  introduced  by  a  low  intensity 

array  of  vortices".  Proceedings  of  the  International  Symposium 
on  Modern  Developments  in  Fluid  Dynamics.  Haifa,  Israel  (16-23 
December  1973),  J.  Rom  (Editor) ,  published  by  SIAM  (1977) ; 
also  SIAM  JL..  APPl-  Math. .  vol.  28,  no.  2  (March  1975) . 

4.  Telionis,  Demetri  P. ,  "Calculations  of  time-dependent  boundary 

layers".  Unsteady  Aerodynamics,  Vol. I,  Proceedings  of  a 
symposium  held  at  the  Univ.  of  Arizona,  edited  by  R.B.  Kinney 
(July  1975) . 

5.  Loehrke,  R.I.,  Morkovin,  M.V. ,  and  Fejer,  A. A.,  "Review  -  Transition 

in  nonreversing  oscillating  boundary  layers",  vol.  97,  no.  4, 
Journal  of  Fluids  Engineering  (December  1975) . 

6.  Rogler,  Harold  L.,  "Freestream  vorticity  disturbances  adjusting  to 

the  presence  of  a  plate  -  -  a  quarter-plane  problem".  Journal 
fi£  Applied  Mechanics,  vol.  44,  no.  4  (December  1977). 

7.  - "The  interaction  between  vortex-array  representations  of 

freestream  turbulence  and  impermeable  bodies",  Paper  75-116, 
13th  AIAA  Aerospace  Sciences  Meeting  (January  1975) . 

8.  —"The  interaction  between  vortex-array  representations  of 

freestream  turbulence  and  semi-infinite  flat  plates1',  Journal 
Of.  Fluid  Mechanics .  vol.  87,  part  3  (1978) . 

9.  - "Rotational  and  irrotational  freestream  disturbances  interacting 

with  a  semi-infinite  plate",  Eighth  U.S.  National  Congress  of 
Applied  Mechanics,  Univ.  of  Cal.  at  Los  Angeles  (26-30  June 
1978)  . 


10.  Rogler,  Harold  L.  and  Reshotko,  Eli,  "Freestream  fluctuations 

interacting  inviscidly  with  semi-infinite  plates" 

(to  be  published) . 

11.  Tollmien,  W. ,  "Uber  die  Entstehung  der  Turbulenz",  Geael 1 schaf fc  der 

Wiasenschaften.  Gottingen.  Mathematisch- 

Natarwissenschaftliche  Klasse.  Nachrichten  (1929). 


97 


AEDC-Tfl-83-4 


12.  Mack,  Leslie  M.r  "A  numerical  study  of  the  temporal  eigenvalue 

spectrum  of  the  Blasius  layer",  jL.  Fluid  Mfifih . ,  vol.73  (1976). 

13.  Corner,  D.,Houston,  D.J.R. ,and  Ross,  M.A.S.,  "Higher  eigenstates  in 

boundary-layer  stability  theory",  i*.  Fluid  Hech. ,  vol.77,  part 
1  (1976). 

14.  Grosch,  C.E.  and  Sal wen,  H.,  "The  continuous  spectrum  of  the 

Orr-Sommerfeld  equation.  Part  1.  The  spectrum  and  the 
eigenfunctions",  i*.  Fluid  Mech. .  vol.  87  (1978) . 

15.  Rogler,  Harold  L.,  "Vorticity  fluctuations  propagating  downstream 

and  interacting  with  Falkner-Skan  boundary  layers" 

(to  be  published) . 

16.  Rogler,  H.  and  Tsuge,  S.,  "Standing  waves  which  decay  or  grow 

exponentially  in  the  streamwise  direction".  Bullet in  of.  the 
American  Physical  Society,  vol. 25,  no. 9  (1980). 

17.  - "Initial-value  problem  and  solutions  for  three-dimensional 

disturbances  in  boundary  layers  and  free-shear  layers". 
Bulletin  of  the  American  Physical  Society,  vol. 26,  no. 9  (1981). 

18.  Rogler,  H.,  "The  production  of  higher  and  lower  modes  and  the  decay 

of  an  array  of  vortices  near  a  boundary".  Bulletin  of.  the 
American  Physical  Society,  vol. 20,  no. 9  (1975). 

19.  Hultgren,  Lennart  S.  and  Gustavsson,  L.  Hakan,  "Algebraic  growth 

of  disturbances  in  a  laminar  boundary  layer”.  Physics  of. 
Fluids .  vol. 24,  no. 6  (June  1981). 

20.  Bechert,  D.  and  Michel,  U.,  "The  control  of  a  thin  free  shear 

layer  with  and  without  a  semi-infinite  plate  by  a  pulsating 
flow  field",  Acustica.  vol. 33,  no. 5  (1975)  . 

21.  Dovgal’,  A.V. ,  Kozlov,  V.V. ,  and  Levchenko,  V.  Ya.,  "Experimental 

investigation  into  the  reaction  of  a  boundary  layer  to  external 
periodic  disturbances" .Fluid  Mechanics  (Russian),  vol. 15,  no. 4 
(July-August  1980);  Translated  from  Izvestiya  Akademii  Nauk 
SSSR,  Mekhanika  Zhidkosti  i  Gaza,  no. 4  (July-August  1980). 

22.  Reshotko,  Eli,  "Remark  on  Engineering  Aspects  of  Transition", 

Transition  and  Turbulence .  edited  by  Richard  E.  Meyer, 
Academic  Press  (1981). 

23.  Nachtsheim,  Phillip  R.  and  Swigert,  Paul,  "Satisfaction  of 

asymptotic  boundary  conditions  in  numerical  solution  of  systems 
of  nonlinear  equations  of  boundary  layer  type",  NASA  E-2959 
(1965) . 


98 


AEDC-TR-83-4 


24.  Prandtl,  L.r  "Bemerkungen  uber  die  Entstehung  der  turbulenz",  Z. 

angew.  Math.  Mach. .  1,  pp. 431-436  (1921) . 

25.  Rockwell,  D. ,  "Oscillations  of  impinging  shear  layers",  AIAA  20th 

Aerospace  Sciences  Meeting,  AIAA  Preprint  82-0047  (11-13 

January  1982) . 

26.  Tsuge,  S.,  "Methods  of  separation  of  variables  in  turbulence 

theory",  NASA  CR-3054  (1978). 

27.  Rogler,  Harold  L. ,  "Waves  which  travel  upstream  in  boundary 

layers",  AEDC-TR-83-4  (to  be  published). 

28.  Ludloff,  H.F.,  De  Santo,  D.F.,  and  Parthasarathy,  R. ,  "Attempts  at 

derivation  of  transition  from  laminar  into  turbulent  flow  along 
a  flat  plate",  AFOSR  TN  60-105,  New  York  University  (January 
1960) . 

29.  Wray,  Alan  and  Hussaini,  M.Y.,  "Numerical  Experiments  in  boundary 

layer  stability",  AIAA  18th  Aerospace  Sciences  Meeting  (14-16 
January  1980) . 

30.  Leonard,  A.,  "Numerical  simulation  of  Interacting, 

three-dimensional  vortex  filaments".  Lecture  Notes  in  Physics , 
vol.  35,  Springer-Verlag,  Berlin  (1975) . 

31.  Grosch,  Chester  E. ,  "Numerical  simulation  of  instability  and 

transition",  Low-Speed  Boundary  Layer  Transition  Workshop 
(16-17  July  1974) ;  Abstract  appears  in  R-1752-ARPA/ONR  by 
William  S.  King  (June  1975). 

32.  Fasel,  H.,  "Investigation  of  the  stability  of  boundary  layers  by  a 

finite  difference  model  of  the  Navier  Stokes  equations", 
Fluid  Mech. .  vol.  78,  part  2  (1976). 

33.  Miller,  G.  and  A.  Callegari,  A.,  "The  effects  of  acoustical 

disturbances  on  boundary  layer  transition",  NYU/DAS  78/01,  New 
York  University/Dept,  of  Applied  Science,  26-36  Stuyvesant 
Street  (Barney  Bldg),  New  York,  New  York  10003  (January  1978). 

34.  Orszag,  S.A. ,  "Turbulence  and  Transition:  A  Progress  Report", 

Proceedings  of  the  Fifth  International  Conference  on  Numerical 
Methods  in  Fluid  Dynamics  Springer-Verlag,  Berlin  (1976) . 

35.  Murdock,  John  W.,  "Tollmien-Schlichting  waves  generated  by  unsteady 

flow  over  parabolic  cylinders",  Paper-AIAA-81-0199,  AIAA  19th 
Aerospace  Sciences  Meeting  (12-15  January  1981). 

36.  Jordinson,  R. ,  "Spectrum  of  eigenvalues  of  the  Orr-Sommerfeld 

equation  for  Blasius  flow".  Physics  of  Fluids,  vol. 14  (1971). 


99 


AEDC-TR-83-4 


37.  Antar,  Basil  N.  and  Benek,  John  A. ,  "Temporal  eigenvalue  spectrum 

of  the  Orr-Sommerf eld  equation  for  the  Blasius  boundary  layer". 
Physics  of  Fluids,  vol.  21,  no.  2  (February  1978) . 

38.  Rogler,  Harold  L.,  "Fluctuations  in  a  boundary  layer  introduced  by 

traveling-wave  irrotational  freestream  disturbances".  Low-speed 
Boundary- Layer  Transition  Workshop:  II,  RAND  Corporation  (13-15 
September  1976) . 

39.  -  "The  coupling  between  freestream  disturbances,  driver 

oscillations,  forced  oscillations,  and  stability  waves  in  a 
spatial  analysis  of  a  boundary  layer”.  Conference  Proceedings 
Of.  the  Fluid  Dynamics  Panel  Symposium  on.  Laminar-Turbulent 
Transition.  Advisory  Group  for  Aerospace  R  and  D  Conference 
Proceedings  No.  224  (2-4  May  1977) . 

40.  Murdock,  John  W,  and  Stewartson,  Keith,  "Spectra  of  tne 

Orr-Soramerf eld  equation",  Physics  fif.  Fluids,  vol. 20,  no. 9 
(September  1977) . 

41.  Salwen,  H.  and  Grosch,  C.  E. ,  "The  spatial  continuum 

eigenfunctions  of  the  Orr-Sommerfeld  equation",  Bull.  Am. 
Physical  Soc..  22,  no.  10  (1977). 

42.  Ellinwood,  J.  W. ,  'Forced  response  of  a  Blasius  flat-plate 

boundary  layer  to  an  external  vortex  street",  SAMSO-TR-79-65  (5 
July  1979) . 


100 


AEDC-TR-83-4 


Appendix  A.  ANALYTICAL  SOLUTIONS  FOR  STANDING  WAVES 
IN  A  UNIFORM  MEAN  FLOW  NEAR  A  WALL 

A.l  INVISCID  SOLUTION 


For  an  inviscid  flow  with  U  =  1,  eqn.(2.1)  reduces  to 


1 


1 

az4 


A  l  v*v  =  o 

dx  j 


and  the  Orr-Sommerfeld  equation,  eqn.(2.6b). 


reduces  to 
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/  — 


+ 
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(A.l) 


(A.  2) 


where  r  =  +1  for  a  growing  wave  and  -1  for  a  decaying  wave.  For  tne 
cases  of  interest  with  p  and  w  both  real  and  nonzero,  the  leading 
coefficient  cannot  vanish.  The  solution  of  the  equation  is  generally 


(A.3) 


where  A  and  B  are  complex  constants.  To  satisfy  the  impermeability 
condition,  then  A--B  and  the  solution  is 

< fztZAesthpy  (a.  4) 

Hence  for  the  inviscid  case  with  a  uniform  mean  flow,  only  a  sine  wave 
arises,  and  the  complex  amplitude  ^is  the  same  for  both  growing  and 
decaying  waves.  This  flow  is  irrotational. 


A. 2  VISCOUS  SOLUTIONS 


Eqn.(2.1)  reduces  to  the  following  equation 

flow 


for  a  uniform  mean 


(A.  5) 


The  Orr-Sommerfeld  equation  reduces  to 


7 


(A. 6) 


The  characteristic  length  5  in  the  nondimensionalizations  is  a  dummy 
parameter  for  this  uniform  flow  case.  While  that  length  could  be 
defined  as  the  characteristic  viscous  length,  the  same 
nondimensionalizations  will  be  used  here  as  defined  in  Chapter  2  so 
that  our  analytical  results  can  be  compared  with  the  numerical 
solutions.  For  solutions  of  exponential  form,  the  four  independent 
solutions  are 
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<p(y)  -  Ae 

where  the  exponent  is 


'/z 


m 


+ jfe  +  ^ cos^y  ¥■  P  f/xjsy 

_ _ _ _  ^ 

-('■p«t-pz-<'**s)'lz = 

Introducing  the  impermeability  and  no-slip  conditions,  then 

<f(y)  =  -  e  ^  +  cos  ^y  -  ~  tt'n^y 

The  choice  of  growing  or  decaying  waves  influences  the  solution  tnrough 
the  exponent  m. 

Rather  than  normalizing  the  solution  by  setting  the  constant 
C  =  1,  the  choice  D  =  1  is  also  valid.  With  this  normalization,  the 

solution  is 


(A. 7) 


(A.  8) 


(A.  9) 


<P(y) 


=  5  -&  cos  8y  +  5<v? 


m 


I  ~ 


As  ,  m  —  eo ,  and  the  inviscid  solution  is  recovered  except  for 

the  arbitrary  constant.  This  normalization  is  not  used  below. 

From  the  continuity  equation,  the  longitudinal  velocity  is 

U  =  /((/)  erf*~'eejt  where  /(y)  =  r  (-  &  e  ^  S/h^y* 

From  the  definition  of  vorticity,  ,  . 

,  rex—  ta/t 

=  ? (y)  e  f  (a. id 

the  vorticity  amplitude  is 

/  (rct>  \  -  mcJ 

Z(y)  =  *s(-/+-p~)e 

Since  m-»<*>  as  Rs— ,  the  rotational  region  is  reduced  to  a  layer  of 
vanishing  thickness  with  infinite  vorticity.  Hence  a  vortex  sheet  is 
positioned  at  the  wall.  For  finite  Reynolds  numbers,  note  tnat  the 
value  of  vorticity  at  the  wall  is  proportional  to  the  Reynolds  number. 

If  the  no-slip  condition,  4L  (0)  -  0,  is  replaced  by  a  no-shear 
condition,  <Puu  (0)  -  0,  then  the  solution  is  ■  A  =  B  =  C  =  0,  or 
^(y)  -  Dsin  y.  Hence,  the  viscous  solution  reduces  to  the  inviscid 
form  for  all  Reynolds  numbers  with  the  no-shear  form  of  the  wall 
boundary  condition.  No  sublayer  or  rotational  region  forms. 
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Appendix  B.  SUBLAYER  SOLUTION  FOR  STANDING  WAVES 


When  Prandtl  (Ref .13)  analyzed  the  viscous  sublayer  and  first 
demonstrated  how  viscosity  could  be  destablizing  in  a  boundary  layer, 
he  made  three  assumptions: 

(1)  The  fluctuating  pressure  is  impressed  upon  the  thin  sublayer. 

(2)  The  mean  velocity  is  sufficiently  small  in  the  sublayer  such  that 
it  can  be  neglected. 


(3)  The  fluctuating  velocity  perpendicular  to  the  plate  can  be 
neglected  in  the  x-momentum  equation. 


These  assumptions  will  now  be  applied  to  the  sublayer  of  the 
standing  wave.  The  objective  is  to  obtain  a  closed-form  solution  which 
illustrates  the  essential  features  of  the  sublayer.  Introducing  these 
assumptions  and  using  the  same  nondimensionalizations  as  defined  in 
Chapter  2,  then  the  x-momentum  equation  reduces  to  the  unsteady,  forced 
diffusion  equation 


^ - 'jf  +  € 


where  if  -  )  6 


rpx- tat? 


(B.l) 


Outside  of  the  sublayer  but  not  too  far  from  the  wall,  the  diffusion 
term  is  neglected,  and  a  balance  exists  between  the  pressure  gradient 
and  the  unsteady  term 


(B.2) 


where  the  pressure  varies  as 

P  =  vt  e 


X-  tc&f 


(B.3) 


where  r=+l  or  -1  for  growing  or  decaying  waves.  Hence  the  complex 
amplitudes  of  the  pressure  and  velocity  are  related  by 


=  rf  %>  (B.4) 

and  this  pressure  is  assumed  to  be  constant  across  the  sublayer.  The 
momentum  equation  then  reduces  to 


-twf  -  -  rpTT"*  efo 

with  boundary  conditions 


(B.5) 


/  =  O  at1  y-O 


a/?d  £  as 


(B.6a,b) 


Seeking  solutions  of  form  exp (my) ,  then  the  homogeneous  solution  is 


(B.7) 
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while  the  particular  integral  is  a  constant  with  the  value 

^  (B.8) 

Hence  the  solution  for  the  amplitude  of  the  longitudinal  velocity  is 

*(?)  -  4  { /  -  *V] }  <B-9> 

This  amplitude  is  plotted  in  Figure  3.4, 
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English 
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D=d/dy 
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NOMENCLATURE 

coefficient  of  the  mth  Chebyshev  polynomial  in  the 
expansion  of  the  qth  derivative  of  & 

c/^/c/f^  where  a*> 

m  ^-o 

constants  in  the  uniform-flow  solution  of  the 
Orr-Sonunerfeld  equation 

complex  phase  speed 

cosine  and  sine  integrals 

complex  constants  used  in  eqn.(4.7) 

ordinary  derivative  in  the  y-direction 

nth  ordinary  derivative  (dV  -0) 

averaged  disturbance  kinetic  energy  per  unit  mass 

complex  constant  used  in  eqn.(4.10) 

complex  amplitude  of  the  longitudinal  velocity 

value  of  f  at  the  sublayer  edge;  value  in  the 
freestream 

exponent;  index  on  Chebyshev  polynomial 
Falkner-Skan  parameter 
number  of  terms  in  the  series 

m-o 

disturbance  pressure 

integer  with  value  +1  for  a  growing  wave 
and  -1  for  a  decaying  wave 

Reynolds  number  based  on  characteristic  thickness  of  the 
boundary  layer 

real  part  of  the  argument  [1 

time 

mth  Chebyshev  polynomial 


105 


AEDC-TR-83-4 


U/V 


U,V 


•>eo 


X 

X(x,t;£) 

y 

y, 

ye 

z=x+iy 

Z(y> 


disturbance  velocities  in  the  x  and  y  directions 

mean  velocities  in  the  x  and  y  directions 

mean  x-velocity  in  the  freestream.  For  Falkner-Skan 
boundary  layers,  the  notation  U^tx)  emphasizes  that 
this  is  the  mean  velocity  at  the  boundary  layer  edge  at 
the  position  x  downstream  of  the  vertex  of  the  wedge 

coordinate  parallel  to  plate  and  in  streamwise  direction 

function  appearing  in  eqn.(1.16)  and  defined  in  (1.17b) 

cordinate  normal  to  the  plate 

phase  shift  used  in  Section  1.3 

y-value  of  the  boundary  layer  "edge" 

complex  coordinate 

complex  amplitude  of  the  fluctuating  vorticity 


Greek  and  Script 


x-waven umber 


y-wavenumber 


Trfrs  is  the  vertex  angle  of  the  wedge  along  which 
the  Falkner-Skan  boundary  layer  develops 


3yz 


gradient  operator 
Laplacian  operator 

characteristic  thickness  of  the  boundary  layer.  For  a 
Falkner-Skan  boundary  layer,  this  length  is  defined  as 

^  _  /  2_  \  Vz 

characteristic  thickness  of  the  viscous  sublayer 


6*  1/R 


inverse  of  Reynolds  number 


7? 


viscosity 

kinematic  viscosity 


f 


vorticity  vector 
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m&'Xr  +  cft' 

r 

P 

60 

Superscripts, 
(  ) 

* 

e 

r,i 

s 


disturbance  vorticity  in  the  z-direction 
amplitude  of  pressure  disturbance 
mass  density 

complex  amplitude  of  the  normal  velocity 
disturbance 

streamf unction 
frequency 


Subscripts,  and  Miscellaneous  Notation 


(a)  ensemble  average  or  suitable  time  average 

(b)  for  periodic  disturbances,  the  time  average  over 


one 


time  period  , 
A-v/to 


=  tfUVl' 


complex  conjugate 


boundary  layer  edge 

real  and  imaginary  parts  of  a  complex  number  or  variable 


sublayer 

partial  derivative  with  respect  to  x 
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